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TECHNICAL NOTE 4276 

AN APPROXIMATE ANALYTICAL METHOD FOR STUDYING ENTRY 
INTO PLANETARY ATMOSPHERES 
By Dean R. Chapman 

SUMMARY 


The pair of motion equations for entry into an exponential planetary 
atmosphere is reduced to a single, ordinary, nonlinear differential equa- 
tion of second order by disregarding two relatively small terms and by 
introducing a certain mathematical transformation. The reduced equation 
includes various terms, certain of which represent the gravity force, the 
centrifugal acceleration, and the lift force. If these particular terms 
are disregarded, the differential equation is linear and yields precisely 
the solution of Allen and Eggers applicable to ballistic entry at rela- 
tively steep angles of descent. If all the other terms in the basic 
equation are disregarded (corresponding to negligible vertical accelera- 
tion and negligible vertical component of drag force), the resulting 
truncated differential equation yields the solution of Sanger for equi- 
librium flight of glide vehicles with relatively large lift-drag ratios. 

A number of solutions for lifting and nonlifting vehicles entering 
at various initial angles also have been obtained from the complete non- 
linear equation. These solutions are universal in the sense that a single 
solution determines the motion and heating of a vehicle of arbitrary 
weight, dimensions, and shape entering an arbitrary planetary atmosphere. 
One solution is required for each lift-drag ratio. These solutions are 
used to study the deceleration, heating rate, and total heat absorbed for 
entry into Venus, Earth, Mars, and Jupiter. From the equations developed 
for heating rates, and from available information on human tolerance 
limits to acceleration stress, approximate conditions for minimizing the 
aerodynamic heating of a trimmed vehicle with constant lift-drag ratio 
are established for several types of manned entry. A brief study is 
included of the process of atmosphere braking for slowing a vehicle from 
near escape velocity to near satellite velocity. 


INTRODUCTION 


One of the many challenging problems connected with space flight 
occurs during the terminal phase of operation when a vehicle at near 
orbital velocity enters the earth’s atmosphere or the atmosphere of 
another planet. Some important aspects of this problem are the possibly 


2 


NACA TN U 276 


severe decelerations for human occupants , the intense aerodynamic heating, 
and the tactical aspect of having satisfactory control over both the time 
and location of landing. The problem is made more interesting by inter- 
relationships between these aspects which require, as always, keen under- 
standing in order to make the best design compromises. For example, the 
lowest heating rates and smallest decelerations are obtained with very 
shallow entry paths, such as would be obtained by letting an orbit 
gradually decay; but the tactical aspects of fixing the time and location 
of a vehicle upon landing are most difficult with these very shallow 
re-entries. Also, the total heat absorbed during descent is greater for 
shallow entries than for steep ones. If descent at a steeper angle is 
induced by deflecting the orbit, such as by means of a retrorocket, then 
the total heat absorbed for laminar flow is reduced substantially, and 
the time and location aspects of recovery are improved, but both the 
deceleration and the heating rate are increased. In order to devise an 
efficient method of entry for a given application, it is highly desirable 
that a designer have available relatively simple equations for computing 
how each variable at his disposal affects the entry trajectory, the 
deceleration, and the aerodynamic heating. 

For several special types of entry, analytical theories are available 
which provide simple equations showing clearly how each variable affects 
the motion and aerodynamic heating. In the case of ballistic-type entry 
without lift at sufficiently steep angles that the gravity and centrifugal 
forces can be disregarded, the analysis of Allen and Eggers (ref. 1) 
provides such equations . In the case of smoothly gliding- type entry at 
zero initial angle with a sufficiently large lift-drag ratio that the 
vertical acceleration and the vertical component of drag force can be 
disregarded, the analysis originally given by Sanger (refs. 2 and 3) would 
be applicable. In the case of skipping vehicles entering at sufficiently 
steep angles and with a sufficiently large lift-drag ratio that the 
gravity and centrifugal forces can be disregarded, the analysis of Eggers, 
Allen, and Neice (ref. 4) would apply. For more general types of entry, 
though, where the gravity force, centrifugal force, lift force, vertical 
acceleration, and vertical component of drag are all of importance, these 
existing analyses would not apply. Such would be the case, for example, 
for the entry of a satellite with a small lift-drag ratio, or for the 
entry of any orbiting vehicle starting with a very small initial angle . 

As a result, present understanding of the relatively shallow entries 
which are of special interest to manned space flight - is based primarily 
on numerical calculations made with computing machines in connection with 
relatively specific vehicles (see, e.g., refs. 5 > 6 , and 7) • 

The objective of the present report is to develop an approximate 
analytical solution to the motion equations which is usable for engineer- 
ing calculations and which is applicable to an arbitrary planetary atmos- 
phere, to a lifting or nonlifting vehicle, and to entries along either 
shallow or steep descents. Such a solution could be applied to a fairly 
broad variety of vehicles, such as skip, glide, satellite, ballistic, 01 
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escape vehicles undergoing the process of atmosphere braking. An addi- 
tional objective is to develop a method applicable to composite types of 
entry, such as entering initially with zero lift, and then suddenly chang- 
ing the lift and/or drag at any number of points during the descent. 

During the preparation of this report an interesting report by Gazley 
(ref. 8) became available in which he considers the entry of a nonlifting 
satellite into a planetary atmosphere from a decaying orbit. He obtains 
an approximate analytic solution by making an arbitrary assumption about 
the relationship between velocity and angle of descent which is not made 
in the present report. As a result, his end equations for this particular 
type of entry are quantitatively different, though qualitatively similar 
to those of the present report, as discussed briefly later. 


NOTATION 


a resultant deceleration 


A reference area for drag and lift, sq ft 

C dimensional constant in heat-transfer equations 

(17,000 Btu ft -3 / 2 sec -1 for numerical calculations of this report) 

Cp drag coefficient, ? — 2 

S Pco^A 


C L 

D 

8 

Sc 

k i 


lift coefficient, = 

!p„v*a 

drag force, lb 

gravitational acceleration, ft sec -2 

gravitational conversion constant, 32.2 ft sec“ 2 

ratio of local heat flux to that at a stagnation point, -3- 


k 2 average value of heat flux relative to stagnation point value , 



l characteristic length of vehicle, ft 
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lift force, lb 

mass of vehicle, slugs 

mean molecular weight of planetary atmosphere (consistent units 
with gas constant and g) 

Prandtl number 

convective heat-transfer rate per unit area, Btu/sq ft sec 
total convective heat absorbed, JJq dt dS, Btu 

dimensionless function proportional to heating rate (u for 
laminar flow) 

dimensionless function proportional to total heat absorbed 
u 3/2 Z~ 1/2 cos~ 2 cp du 


distance from planet center, ft 

universal gas constant, or radius of curvature of vehicle surface 
in feet 


Pfl VI 

Reynolds number, 

M-oo 

circumferential distance traveled, ft 
surface area wetted by boundary layer, sq ft 
time, sec 

temperature (various units employed) 

circumf erential velocity component normal to radius vector, ft/sec 
circular orbital velocity, ■Jgr, ft/ sec 

ratio, — 
u c 

upper limit for range and total heat absorbed (see eqs . (28) and 
(39b)) 


altitude, ft 


NACA TN ^276 


5 


v vertical velocity component (along direction of radius vector), ft/sec 
V resultant velocity, n/u 2 + v 2 

W weight of vehicle at earth T s surface, mg c , lb 

Z dimensionless function of u determined by equation (2l) and 
appropriate boundary conditions 

3 atmospheric density decay parameter, ft " 1 

7 ratio of specific heats behind bow wave 

6 angle in polar coordinates 

[x coefficient of viscosity, slug ft" 1 sec " 1 

p density, slug ft " 3 

cp flight-path angle relative to local horizontal direction; positive 
for climbing flight, negative for descent 

Subscripts 

o sea level 

00 free stream 

s stagnation point 

1 initial condition 

W 

b break where is dis continuously changed 

e relative to earth 

Superscripts 

r differentiation with respect to u 

mean value for exponential approximation to atmosphere density- 
altitude relationship, or dimensionless quantity 
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, ANALYSIS 

Assumptions and Approximations 


The problem analyzed concerns that portion of the descent of a 
vehicle into a planetary atmosphere wherein the decelerations and the 
convective aerodynamic heating are dominant. Three assumptions made at 
the outset are: 

(i) Atmosphere and planet are spherically symmetric. 

(ii) Atmosphere density p ro varies exponentially with altitude. 

(iii) Peripheral velocity of planet is negligible compared to the 
velocity of the entering vehicle. 

Assumption (i) is reasonable for those planets which have only small 
equatorial bulges (such as Venus, Earth, and Mars), inasmuch as the severe 
aerodynamic heating and decelerations occur over a length of flight path 
which is small compared to the planet's mean radius (the order of 
one tenth the planet radius for nonlifting bodies) . The assumption of 
spherical symmetry, however, would not be as reasonable for planets with 
relatively large equatorial bulges, such as Jupiter and Saturn. As noted 
later, this assumption of spherical symmetry can introduce some inaccuracy 
if the descent is nearly along a line of longitude and if the vehicle also 
happens to have a relatively large lift-drag ratio . For large lift-drag 
ratios the important deceleration and heating portions of the descent can 
be prolonged over a distance comparable to the planet's radius; hence, 
the nonspherical nature of the atmosphere could be important in such 
cases . 

Assumption (ii), of an exponential atmosphere, is based upon the 
simple kinetic theory of an isothermal gas in a uniform gravitational 
field. This theory yields the well-known exponential approximation for 
atmospheres (see ref. 9 , ch. Ill, for example) 



( 1 ) 


where 



( 2 ) 


and where M is the mean molecular weight of the planet's atmosphere, 

T the mean temperature, R the universal gas constant^ and g the local 
acceleration due to gravity. It is to be noted that 'p 0 represents the 
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intercept of the straight line which best fits a curve of log p versus 
altitude, and is not the same as the true ‘sea level density p Q . From 
data such as presented in references 10 and 11, approximate mean values 
of several quantities of interest for various planets are as follows 
(the subscript & designates a value relative to the earth): 


Planet 

r a> 

g® 

Gases 

Pe 

M, 

gm mol -1 

T, 

°K 

P"S ft 

■v/(0r)* 

Venus 

0.97 

0.87 

co 2 ,n 2 

0.8 

4o 

270 

2x10 4 

1.0 

Earth 

1.0 

1.0 

n 2 ,o 2 

1.0 

29 

240 

2.35X10 4 

1.0 

Mars 

.53 

.38 

n 2 ,co 2 

1.0 

30 

220 

6x10 4 

.47 

Jupiter 

11.0 

2.63 

h 2 ,ch 4 

.5 

3 

170 

6x10 4 

2.0 


The exponential approximation for the earth (with 'p 0 = 0.0027 slug ft -3 ) 
is compared in figure 1 with the relatively recent (1956) ARDC model of 
the atmosphere. It is evident that a single value for 3 appears to be 
a reasonable approximation at altitudes below about 400 , 000 feet (80 miles, 
roughly) . In most cases peak decelerations and maximum aerodynamic heating 
occur well below this altitude. Moreover, the region of most important 
heating and deceleration for a given vehicle occurs only over a relatively 
thin strip of altitude (very roughly over a 70,000-foot strip across which 
the density changes by about a factor of 20) . Since the analysis which 
follows enables the altitude of this important strip to be calculated 
quickly for any given vehicle, the exponential decay parameter 3 in each 
case could be selected, if desired, as corresponding to this particular 
altitude rather than to the mean value tabulated above . A plot of the 
dimensionless parameter 3r as a function of altitude is shown in fig- 
ure 2 for the ARDC model atmosphere. In determining 3r consideration 
is given only to the 70,000-foot region of air immediately above a given 
altitude. The fluctuations in \f$r for this standard atmosphere below 
about 400,000 feet amount to the order of ±10 percent from a mean value 
of 30 and are attributed primarily to the variation in temperature with 
altitude . Inasmuch as variations in temperature with season and with 
latitude (see ref. 12, for example) can fluctuate the order of ±15 percent, 

__ _ 2 _ j 2 

the parameter T can fluctuate about +7 percent. For most numer- 

ical calculations in this report, a constant value «/]3r = 30 is used for 
the earth's atmosphere corresponding to a mean atmospheric temperature of 
240° K (432° R) . 

Assumption (iii), that the peripheral velocity of the planet is 
negligible compared to the velocity of the entering vehicle, would not 
introduce significant errors for most descents into most planetary atmos- 
pheres . For descents nearly along a line of longitude, the errors in heat 
transfer and deceleration would, of course, be negligible. The greatest 
error would occur in an equatorial descent. As a measure of this error, 
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we can take the ratio of the equatorial peripheral velocity u p of the 
planet to the circular satellite velocity u c . This ratio for several 
planets is as follows: 


Hence the error introduced by assuming a nonrotating atmosphere in the 
case of near-equatorial descents would be negligible for Venus, appreci- 
able though not large for Earth and Mars, but probably significant for 
Jupiter . 

In addition to these three physical assumptions, two mathematical 
approximations are made in the development of the subsequent analysis in 
order to effect major simplifications in the structure of the equations 
of motion. They are mentioned here for convenience: 

(a) In a given increment of time, the fractional change in distance 
from the planet center, dr/r, is small compared to the fractional 
change in velocity du/u; that is, |dr/r|<< |du/u|. 

(b) For lifting vehicles, the flight-path angle cp relative to the 
local horizontal direction is sufficiently small that the component 
of lift in the horizontal direction is small compared to the drag; 
that is, |(L/D)tan cp ) << 1. 

For nonlifting vehicles (e.g., ballistic entry), approximation (b) 
is automatically satisfied; approximation (a) does not specifically 
restrict the descent angle (0 < “ > to 90 1 ” 1 can be analyzed for nonlifting 
vehicles), but it does restrict the analysis to a portion of the over-all 
trajectory below an upper altitude limit. Above some altitude dr/r 
cannot be small compared to du/u, as is shown to be the case on mathe- 
matical grounds in appendix A. Physically, this is clear from the law 
of conservation of angular momentum which states that in the absence of 
drag, d(mur) = 0, or dr/r = -du/u. Consequently, the present solution 
would be reasonable at least below an altitude where drag has slowed down 
a vehicle slightly to some point (A in sketch) where dr/r = 0.1 du/u. 

It is shown in appendix B that this corresponds to tne point where drag 
has reduced the vehicle velocity by about 0.01 of the initial velocity. 
Above this altitude (point A), orbit-type calculations could be applied. 

A method for joining the present solution to Keplerian ellipses is 

discussed in appendix B. 

1 This value may be a factor of ten higher due to the uncertainty in 
the length of the Venus day. 



Venus 

Earth 

Mars 


Jupiter 


x 0.002 

.06 

.07 

.29 
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For lifting vehicles (e.g., skip or gliding entry) assumption (b) 
clearly restricts the analysis to small angles of descent . Even if a 
lifting vehicle starts entry horizontally, the angle of descent will 
increase as the velocity is reduced (and as the centrifugal forces are 
diminished) until (L/D)tan cp becomes unity in the terminal subsonic 
gliding phase. Although the solution is not valid, strictly speaking, 
when (L/D)tan|cp| is comparable to unity, a reasonable over-all trajectory 
would be obtained by stopping the present solution at the point where 
-(L/D)tan 9=1 (point B in above sketch), and considering that 
-(T./p) t an cp = 1 thereafter. As sketched, peak heating and maximum decel- 
eration occur well within the range (solid line) where the present solu- 
tion applies . 

The limitations resulting from approximations (a) and (b) are 
examined in appendix A, where it is shown that for vehicles entering from 
decaying satellite orbits, with or without positive lift, the errors 
introduced are only the order of a few percent insofar as aerodynamic 
heating and peak decelerations are concerned. Surprisingly small errors 
result from approximation (b), even for very large l/d ratios, because, 
in orbital decay or in a smooth glide , the larger the L/D the smaller 
the angle cp at conditions near maximum heating and peak deceleration; 
this keeps the product (L/D)tan 9 small. 

Various modes of entry and the portions of the trajectories of 
satellite, ballistic, escape, glide, and skip vehicles to which the anal- 
ysis applies are sketched in figure 3* 


10 


NACA TN 1+276 


Development of Differential Equation 


Descent in a spherically symmetric atmosphere about a spherically 
symmetric planet would occur in a meridian plane in the absence of lateral 
forces . This confines the problem to one of two dimensions for which 
polar coordinates (r, 0) are convenient. The velocity components are 
(v, u) , respectively, as sketched below. 



The vector acceleration in terms of the unit vectors e r and eg for polar 
coordinates is 




( 3 ) 


where e r and eg are the unit vectors in the r and 0 directions, 
respectively. The local flight-path angle cp (negative for descent) is 


tan cp = — 
u 


(*0 


The vector aerodynamic force 
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must equal the mass m times the vector acceleration in the absence of 
thrust-type forces. Hence, equations (3)> (4) , and (5) yield two com- 
ponent equations of motion 


- = - — = g - — - ^ cos cp + 5 sin cp (6) 

dt 2 dt r m m 


du 1 uv __ D 
dt r m 


cos 


m ^ L 
Cp+ D 


sin cp 


(7) 


It is noted that g and r are local values in these equations. 

We will solve this system of equations by disregarding the term 
uv/r in equation (7) (which, as will be evident shortly, is equivalent 
to assumption (a) that | dr/r | << | du/uj ) . This restricts the solutions 
to problems wherein |uv/r | « | du/dt | , but the restriction is not serious 
for the aerodynamic heating and deceleration aspects of entry. In the 
case of orbital entry, for example, maximum deceleration and heating 
occur at such small angles that uv/r is the order of 1 percent of du/dt 
(see appendix A). An alternate view of what the approximation involves 
can be seen as follows: 



« 


1 


( 8 ) 


Consequently, the disregard of uv/r is precisely equivalent to approx- 
imation (a) mentioned earlier; namely, that the percentage change in 
distance from the planet center is small compared to the percentage change 
in velocity. We will employ this approximation several times more in the 
analysis. Inasmuch as du/u is relatively large only when the drag is 
important, it is understandable why the basic approximation | dr/r |« |du/u| 
yields results applicable to regions of important deceleration and aero- 
dynamic heating, but not to the outer regions of space where orbit-type 
calculations (which do not neglect the acceleration term uv/r compared 
to du/dt) are necessary to describe the motion of a vehicle. In these 
outer regions, radiant heat dominates, while convective heating and 
deceleration are very small. 

By utilizing approximation (a) (inequality (8)), we have 

= - ^ cos cp ( 1 + ^ tan cp 


du 

dt 


(9) 
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so that, by introducing the drag coefficient, the exponential approximation 
to the atmosphere, and approximation (b) (|(L/D)tan cp | « l) , and noting 
that V = u/cos cp, 


— -By 

dU = _ Po e U 2 

dt / m\ cos cp 

W 


We will select as an independent variable 

u = — = u 
u c */gr 


( 10 ) 


( 11 ) 


representing the ratio of horizontal velocity to the local circular 
satellite velocity. The basic approximation (8) taken together with the 
relation dg/g = -2dr/r resulting from Newton's gravitational law enables 
us to disregard derivatives of both g and r relative to derivatives of 
either u, or u; for example, 

du = d( sfgr u) = j— au ( 12 ) 

dt ~ dt dt 


By introduction of the drag coefficient, the motion equation (6) for an 
exponential atmosphere becomes 


1 dv 
g dt 


1 dfy _ -2 Po °D Aru2 e" |3y 
g dt 2 U + 2 m cos 2 cp 


sin cp 


L 

— cos 
D 



(13) 


In order to reduce the pair of motion equations (10) and (13) to a 
single equation, we transform to a new dimensionless dependent variable 
Z defined by 


Z = 



(14) 


and employ u as the independent variable. 2 Thus, by differentiating 
(Z' = dz/du) and keeping in mind the basic approximation (8) 

2 The author knows of no simple way to explain a priori why this coor- 
dinate system Z(u) should be introduced. It was discovered by trial and 
error after trying various other transformed coordinate systems which did 
not reduce the pair of motion equations to a single equation. 
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Z* 


fit 


U U 2 


2(12-' 

W. 


>-[3y 

du 


_ Z dy dt 
u dt du 


(15) 


We see from equations (10) and (12) that 


22= -TiP 

dt cos cp 


uZ 


(16) 


so that substitution into equation (15) and noting that 
dy/dt = v = u *J~gr tan 9 yields 


. Z = /I eo^jp = ^ 

u vg u dt . 


(17) 


Proceeding now by differentiation of v and sin 9 from equation (17)> 
there results 


1 dv _ R d 
e dt “ v g dt 


g dt v S dt \ cos 9 J dt \cos 9 


u_sin_ 9 \ _ 1 du 


1 du f uZ" , u -/ftr sin 2 9 dqA 

du^ 


cos 2 cp 


(18) 


The term d 9 /du representing flight-path curvature can be expressed in 
several ways in terms of the Z function by noting from equations (17) 
and ( 12 ) that 


u -4 (z 1 - £) = u^i. s ^5L 

du \ u/ du 


= uZ" - Z» + = 
u 


= uZ" - ’Tpr sin 


Alternate forms of 
} terms representing (19) 
flight -path curvature 
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Consequently, we can substitute the first form of this equation, together 
with equation (l6) into equation (18) to obtain 
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1 dv 
g dt 


i 5fZ , _1A_ luZ" + tan 2 * fu i (v 
S dt 2 cos 2 cp l L du \ 



( 18 a) 


We note from equations (l4) and (17) that equation (13) can be written 
in the form 


1 dv 
g dt 


iSSui-n 2 + -55-fz' - I 

s dt 2 cos 2 cp \ u 


JJr ^ cos cp^ 


( 20 ) 


Hence, by comparing this equation with equation (l8a), and by observing 
from the second form of equation (19) that 


Z* 




the final equation for the Z function is obtained. 


u 


_d 

du 




1-u 2 


_ U cos 4 cp + «7 Pr — cos 3 cp = 0 
uZ E 


( 21 ) 


In this equation, cos cp = - sin 2 cp can be expressed in terms of Z 

and Z through equation (17) 


\fpr sin cp = Z* - ^ 


(17) 


Thus, the pair of motion equations has been reduced to a single, second- 
order differential equation by using u as the independent variable and 
Z as the dependent variable. 3 For nonlifting vehicles (l/D = 0) the 
equation is applicable to large angles of descent as well as small . For 
lifting vehicles it is applicable for |(L/D)tan cp|«l. In all cases it 
is applicable when |dr/r |/|du/u|«l. We note from equations (4) and (lb) 
that 


3 Clearly, the same reduction would be achieved by using g(u) as the 
independent variable and Zh(u) as the dependent variable, where g(u) 
and h(u) are arbitrary functions. 


NACA TN 427 6 


15 


dr/r uv/r -u sin cp 

du/u du/dt ^ z ^ 22 ' ) 

As noted in appendix A, the ratio | dr/r | / |du/u| is less than 0.1 below 
the altitude where drag has reduced the velocity by about 1 percent of the 
initial velocity. 

The nonlinearity of equation (21) is due to the term (l - u 2 )cos 4 T/uZ 
which represents the effects of gravity and centrifugal forces in inducing 
a curved flight path. It is noted that the basic equation is independent 
of the physical characteristics Cp, W, A of the vehicle as well as 
independent of the sea-level characteristics and g Q . Aerodynamic lift 
occurs only in the combined parameter L/D. The equation has a singu- 

larity at Z = 0 which must be handled analytically in numerical methods . 

A method of solving this equation is discussed in appendix D. 

It is instructive to consider the physical meaning of each of the 
terms in the differential equation (21) . Equations ( 19 ) and (20) help 
in this regard. 


uZ” 



1-u 2 

uZ 


cos 4 cp 




D 


cos 3 cp 


(21a) 


vertical vertical component 
acceleration of dra g force 
( \ / pr sin cp) 


gravity minus lift force 
centrifugal force 


By understanding the physical significance of the various terms one can 
judge, for example, what terns to consider in obtaining special approxi- 
mate solutions . 


Since the basic differential equation is of second order, we need 
two initial conditions to complete the system. We take these at some 
initial velocity uj, and write as generalized initial conditions 

Z(ui) s Z ± Z*(u i )=Z i t 

If the vehicle starts at a very high altitude where the density is negli- 
gible compared to that near peak heating, then the definition (l4) 


f&c C D A [r\ 
\ 2W v 3 / 


u iPi 


(24a) 


l6 
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shows that Z± is very small in such cases . For simplicity we take 
Z± = 0 for entries starting at very high altitudes. The equation 


Zj_’ = \Tp r sin cp^ + ~ 


u- 


(24b) 


shows that Z± * would he equal to *7 0r sin cp when Z± — 0. As an 
example , entry from a decaying satellite orbit (cpi = 0 and ui = 1 in the 
stage of decay before appreciable aerodynamic heating begins), would be 
represented by the initial conditions 

Zi(l) = 0 Zi’(l) = 0 (25) 


One universal Z function would be required for each value of the 
parameter -JJr L/D appearing in the differential equation (21) . 

By allowing Z± * to take on values other than zero and allowing ui 
to be either less than or greater than unity, we can obtain the corre- 
sponding Z functions for ballistic, glide, skip, or escape vehicles 
entering a planetary atmosphere from very high altitudes. By further 
allowing Z± to be other than zero, corresponding Z functions can be 
obtained for entry starting from an initial altitude where the density 
may not be negligible compared to that near peak heating. Before present- 
ing some solutions to equation (21) , though, it is advantageous to show 
how the Z functions, once computed, can rapidly be used to determine a 
number of useful quantities in practical calculations . 


Summary of Some Useful Quantities Related to the Z Function 


From the Z functions, it is a relatively simple matter to obtain, 
for example, the horizontal component of deceleration a 0 by using 
equations (3)> (12), and (l6) , 

_ au = i-m uz ( 26 ) 

y dt cos cp 
or 

_ — s 30 uZ for Earth, cp small 
g dt 

Strictly speaking, g and r are local values in the outer layers of the 
atmosphere where the deceleration takes place. For Earth, however, these 
are not significantly different from their respective surface-level 
values. Local and surface values might be greatly different, though, for 
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planets such as Jupiter and Saturn which are believed to have a very deep 
atmosphere. The equation for the angle of descent (cp<0 for descent) is, 
from equation (17), 


sin cp = 



Z 1 - (z/u) 

30 


for Earth 


(27) 


The circumferential distance traveled between a point where the dimension- 
less velocity is ui and a point where it is us can be expressed in 
terms of Z from equation (l6) , 


As 

r 



1 r Ul cos cp du 

r J u 2 z 


( 28 ) 


or 



for Earth, cp small 


(29) 


Inasmuch as the analysis is not valid in a very small neighborhood of 
Z = 0 where u = uj[, but becomes valid after drag has reduced u by less 
than 1 percent (as shown in appendix B) , we select an upper limit such as 
ui = 0.995 ui or u x = 0.99 ui for the entry range. In a practical appli- 
cation, this range would have to be joined to the range of the appropriate 
Keplerian ellipse in order to obtain the total range . The corresponding 
time elapsed is obtained also with the aid of equation (l6) 




ds 

u 


'Ul 


'/Pi ^ U; 


u 2 


cos cp du 
uZ 


s 27.0 



sec for Earth, cp small, g = g Q 


(30) 


Another useful quantity is the density ratio, referred to the true sea- 
level density (p 0 = 0.00238 slug ft -3 ), which comes from the defini- 
tion (l4) for Z 
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_ _2_ /p / m \ Z 

Po p o v r W/ u 


(31a) 


or 


30 P, 


s/p r p c 


= 3.7 


W 

c d a 


5 X 10- 5 
u 


for Earth, 


W 

c d a 


in lb ft 


-2 


The left side of this last equation is a function only of the altitude 
for a given atmosphere for the ARDC model atmosphere is shown in 

fig. 2), so that it provides for a given Z(u) function the altitude- 
velocity relationship for any model atmosphere. The density ratio 
referred to the effective sea-level density (Pq = 0.0027) which best fits 
the p(y) curve is 



(31t>) 


= 3.2 



Z 

u 


X 10" 5 


for Earth, 


W 

CpA 


in lb ft -2 , */pr = 30 


The dynamic pressure is 


2 P °° CpA cos2cp 


(32) 


= 10 uZ ( for Earth, — — in lb ft -2 , cp small 

VCd V c d a 


and the free-stream Reynolds number per unit length is proportional to 
Z (eq. 14)) 


R^oo _ V Poo 2 «/g3 


1 P'00 [i cos cp VPjjA 


m 


( 33 ) 


= 7100 ( Z for earth, cp small, — — in lb ft -2 


CnA. 


C d A 
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The viscosity of air at the mean atmospheric temperature T = 432° R is 
employed to obtain the constant in this last equation which is valid only 
for the earth* s atmosphere. 

It is interesting that the Reynolds numbers involved during entry 
from a decaying satellite orbit are relatively small. Near peak heating, 
for example, we will see subsequently that the value of Z ranges from 
about 0.17 to 0.015, for L/D ratios between 0 and 1, so the correspond- 
ing Reynolds numbers are of the order of 1000(w/Cj)A) to 100(w/CpA) per 
foot. These are sufficiently small for one to be optimistic about the 
practical possibilities of maintaining laminar flow for shallow entry 
from a satellite orbit. For steep entries, as for ballistic vehicles, 
the Z function is larger, and hence the corresponding Reynolds numbers 
are larger. Curves illustrating this are presented later. 


Fairly simple expressions also can be obtained for the aerodynamic 
heating rate per unit area (q) and the total heat absorbed per unit area 
Q/S . Following the analysis of Lees (ref. 13) , we will consider the 
heating rate at any point on a body to be a certain fraction 


k i s cf- (3*0 

^•s 

of the heating rate q s at a stagnation point of radius of curvature R. 
The heating rate in hypersonic flow at a stagnation point , can be 
expressed as 


^s 




Btu ft -2 sec -1 


(35) 


where the constants C, n, and m depend on the type of boundary-layer 
flow. For laminar flow we have n = l/2 and from the several references 
listed (with p G being the true sea-level density) 


Reference C m Remarks 

14 16,800 3-1 Intermediate enthalpy theory 

13,7 19,800 3.22 Theory of Lees 

15 17,600 3-15 Correlation of AVCO shock- 

tube experimental results 

We will base all our numerical calculations on laminar flow (n = l/2) , 
and will use the value m = 3 for purposes of simplicity (this corre- 
sponds to a gas with viscosity proportional to T 1 / 2 ) , and the value 
C = 17>000 Btu ft _3/2 sec -1 which is adjusted to match a mean of the 
above results for air at velocities near peak heating (u = 0.8). For 
gases other than air we use the theory of Lees (ref. 13) to obtain for 
hypersonic flow C ~ ■7P 0 l- t 0 u c 3 Pr -2 / 3 [ (7 - l)//] 1 / 4 * In subsequent calcu- 
lations, differences in the Prandtl number and in the ratio of specific 
heats for various planets are disregarded. 
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Proceeding from equations (35) , (3*0 > and (31) with n = l/2, m = 3> 
and C = 17,000, the laminar convective heat-transfer rate can be written 
in terms of the Z function and the relative planetary constants 

(g® = s/g e arth> = ^ dearth' e^ 0 *) as 


q = k 1 q s = 590 jpr“2/3^ o l/2 g 3/2 r 5/4pl/4| 


m 


'CpAR 


cos 3 cp 


Btu 

ft 2 sec 

(36) 


where q = u 5/2 Z l/2 . (if the flow were turbulent q ~ u 2 Z , approxi- 
mately, and we would have different powers appearing in equation (36).) 

It is noted that the different variables affect the heat flux in a form 
represented by a series of factors; the expression in curly braces rep- 
resents the effect on heat flux of the particular planetary atmosphere, 
the expression in brackets represents the effect of the physical charac- 
teristics of the vehicle, that is, the mass, dimensions and shape of the 
vehicle, and the dimensionless function q = d 5/2 Z 1/2 represents the 
effect of the particular type of trajectory as determined by the lift- 


drag ratio. 


Whereas equation (36) for heating rate would be useful in studying 
vehicles designed to operate at radiation equilibrium temperatures, an 
equation for the total heat absorbed during entry is of more interest for 
heat-sink type vehicles. 


where 



(37) 


(38) 


is the factor which takes into account the variations in heat flux over 
the whole surface S wetted by the boundary layer. (For a hemisphere, 
for example, ks = 0.5.) Combining equations (37) > (36), and (30) yields 
the following equation for the heat absorbed between ui and u 


where 


fPr" 2/3 M-o 1/2 gr 5 ' 4 1 

V 

rn 

1 3 1/4 X 


J CpAR 


Q (39a) 


Q 


u i —3/ 2 , — 

u du 


-x^ 


, ^ /2 cos 2 cp 


(39b) 


NACA TN 4276 


21 


Heat radiation from the surface is not considered in these equations. 

They are useful in studying vehicles incorporating heat sinks or ablation 
cooling under conditions when the heat radiated away is small compared to 
the heat absorbed; Q = Zffl'(cm) sink (where c is the effective specific 
heat of the sink material) is proportional to the heat-sink weight. We 
note here that the particular planetary atmosphere (g, r, 0) and espe- 
cially the particular trajectory Z(u, L/D) affect the heating rate q 
in a different fashion than the total heat absorbed Q. Examples 
illustrating this are presented later. 


Some Approximate Analytical Z Functions 
Obtained From Truncated Basic Equation 


By disregarding three different combinations of terms in the basic 
differential equation (21), three special solutions are obtained which 
yield results identical to previous approximate solutions . The details 
are described in appendix C and lead to the following approximate 
solutions : 

Terms disregarded 




Solution 


Vehicle 

(see eq. (21a)) 


Z I 

= «/0r (sin cp)u Zn 

u 

Ballistic 

Gravity, centrifugal 
and lift forces; 
cp = cp = constant 

(40) 

Z II 


1-u 2 



Vertical acceleration 

(41) 

u 

^.(5) 


Glide 

and vertical component 
of drag force; cos cp = 1 

z in 

= u 

r z i 1— 

— +v0r cp,- Zn 

L u i 

u 

Ui 

Skip 

Gravity and centrifugal 
forces; cos cp = 1 

(42) 


ra tl 

(whrj) 

c — » 

1 I 






The Zj function provides an approximate solution for the motion and 
heating identical to the solution of Allen and Eggers (ref. l) for balli- 
stic entry. The Zji function corresponds to equilibrium gliding flight 
originally discussed by Sanger (ref. 2). The corresponding aerodynamic 
heating problems for this type of hypersonic flight have been discussed 
by Eggers, Allen, and Neice (ref. 4) who also obtained a solution equiva- 
lent to the Zjji function for skip vehicles. As will be apparent later, 
the Zj function for ballistic vehicles is quite accurate for angles of 
descent greater than several degrees C/0r|cpi|>2 approximately) and the 
Zjj function for hypersonic glide vehicles is quite accurate for l/d 
ratios greater than about 1 Q fir l/d> 30 approximately) provided cpp = 0. 
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The accuracy of the Zjjj function for skip vehicles, however, depends 
on both L/D and the initial angle cpi- The conditions for applicability 
can be determined from an approximate solution which considers both the 
gravity and centrifugal forces that were neglected in obtaining Zm. 

In appendix C the following approximate solution is developed for satel- 
lite entry (up = 1, Zj. = 0) at small initial angles cpi: 

Zjy = u ./pr -Icpj^Zn u - 


By comparing with Zj j j , we see that the gravity and centrifugal forces 
can be disregarded provided 20r|cpi(L/D)| » 1 . An interesting feature 
deduced from Zjy in appendix C is that the total heat absorbed in the 
first skip (which is perhaps the most important) is essentially independent 
of both the initial angle q>j_ and the velocity of exit from the skip. 

The heat absorbed varies as l/s/Cp and hence is a minimum for entry at 
(see appendix C) . For flat plates in Newtonian flow this corre- 
sponds to an optimum L/D of 0.7* 


1 + 


2prcPi(L/D)_ 


In^u 


(43) 


Some Z Functions Obtained From Full Equation 


Entry from a decaying o rbit for various L/D (up = 1, <Pi = 0) -~ We 
turn now from the special solutions obtained by truncating the lull equa- 
tion (21), to some solutions of the complete nonlinear equation applicable 
to vehicles entering from a decaying satellite orbit. As the apogee of an 
elliptical orbit is slowly reduced by drag (primarily exerted near the 
perigee) , the orbit eventually becomes a near circle and then begins a 
gradually decaying spiral; hence, the initial angle cpi for this type of 
entry is taken as zero, and the initial velocity u* = 1. The peak heat- 
ing and the maximum deceleration occur at such small angles that 
cos cp = 1. The differential equation (21) is then 


u i T— - -V 1 + -ffc £ - 0 (wo 

du \du u/ uZ D 

and the corresponding boundary conditions for decaying orbits are 

z(l) = 0 Z»(l) = 0 (45) 

This system need be solved only once for each value of the parameter 
\Tpr(L/D), and the results are then applicable to any planet and to any 
vehicle with arbitrary shape, size, or mass. In particular, the universal 
Z function for L/D = 0 is presented in figure 4(a) (uZ is plotted since 
this product stays within smaller bounds than Z) . Solutions of equa- 
tion (44) also have been carried out for various values of */0r(L/D). The 
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numerical method employed is described in appendix D. Curves of the 
Z functions and related quantities are presented in figure 4(b) for 
values of n/]3t(L/D) corresponding for earth to L/D = 0 . 1 , 0 . 25 , 0 . 5 / 
and 1. They are plotted in coordinates especially suited for comparison 
with the Zjj function of equation (4l) representing Sanger f s concept 
of an equilibrium hypersonic glide. Values are not shown for L/D > 1 
since Zjj can be used in the velocity range of interest for these cases. 
This is evident from the various curves in figure 4(b) . The dashed curve 
represents Zjl which is the exact solution for L/D = 00. From the Z 
functions the various quantities of engineering interest, such as the 
deceleration, descent angle, range, time, density- velocity relationship, 
dynamic pressure, Reynolds number, heating rate, and total heat absorbed 
can be computed from equations (26) to (39) presented earlier. 

Nonlifting entry with initial angle of descent (ui = 1, 9i < 0).- 

We now consider entry when the initial descent angle is not negligible , 
as it is in the case of a decaying orbit, but is some finite value cpq • 
Entry with an initial angle occurs in the case of a ballistic vehicle, 
or a satellite to which a retrorocket has been applied to divert the orbit 
into one which will induce the entry process. The differential equa- 
tion (21) for nonlifting bodies is applicable for large as well as small 
angle s of de s c ent . 


uife 
du \du 




cos 4 cp = 0 


(46) 


The initial conditions are 

Z(ui) = 0 Z*(ui) = »/3r sin cpi 


(^ 7 ) 


In this case we have a double parameter family of solutions (uj^ and 
Jpr sin cp ). Actually, we need solutions to the nonlinear equation (46) 
only for quite small initial angles inasmuch as the Allen-Eggers solution 
(eq. (40)) is applicable for moderate and large angles. This may be seen 
from figure 5 which presents example Z functions corresponding to the 
nonlinear equation for various -cpj_ up to 20° with u^ = 0.9 (23,400 fps, 
for earth). Since the ordinate is Z//j3r( -sin cp ± the^Allen-Eggers 
solution is represented by the ordinate function u Zn(u^/u) on this plot. 
It is evident that their solution, which neglects gravity and centrifugal 
forces, is quite accurate near peak heating (u = 0.7) for descent angles 
greater than about 5 C . Near maximum deceleration (u = 0.4) the descent 
angle has to be somewhat larger for comparable accuracy. It is clear 
that, as far as peak heating and maximum deceleration are concerned, a 
family of solutions to the nonlinear equation need only be computed for 
small initial angles. 

The Z functions for small initial angles and for the case of 
satellite entry (uq = l) are of special practical interest. These are 
presented in figure 6(a) for various values of J (3r cp-^ such that in the 
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earth's atmosphere -cp-j_ = 1° , 2°, 3°.> 4°, and 6°. Rather than to plot 
Z itself, the quantity 30 uZ is plotted which represents for the earth 
the horizontal deceleration in g's. Tabulated values are presented in 
table I for -cp ± = 0°, 0.5°, 1°, 2 °, 3°, and 4°. It is noted that these 
values tabulated are solutions to equation (46) with the cos^cp term 
included, and hence are applicable to terminal conditions of small u 
(say less than 0.1) where cp is large as well as to conditions near peak 
heating and maximum deceleration where cp is small. The tables of Z 
apply to any planet for the same initial value of (3r cpi • The supple- 
mentary tables of -<Pearth> (As/r) earth > an< ^ ^earth can a PPli ed - ' to 

other planets over the range where cp is small by regarding the tabu- 
lated values as representing -(«Vpr /30)9> / 30) (As/r) , and 27 \[$g t, 

respectively (see eqs. (2j), (28), and ( 30 )). 

Entry with initial angle of descent for various L/d (uj = 1, cpj<0) -~ 
If we now consider a vehicle with lift, we must restrict our considerations 
to initial angles of descent -cpi and to the portion of trajectory 

over which -cp remains sufficiently small that (L/D) |tan cp|«l (assump- 
tion (b)). The basic differential equation (21), with cos cp = 1, becomes 
the same as equation (44) , and the initial conditions are now 

z(l) = 0 Z»(l) = J&r cpi (^8) 

Solutions to equation (44) with these initial conditions have been obtained 
for various values of *] pr cpi and for various values of the parameter 

./F(l/d) . 

In figure 6 some curves representing Z functions are presented 
as a function of -cpi for earth. The various portions of this figure 
correspond to L/D for earth of 0.25, 0.5> 0.7> 3J^d 1. It is evident 
from these figures, as might be expected, that small values of L/D and 
_cp ± do not result in any significant skipping, but once the L/D is 
increased beyond a certain amount, or the initial descent angle is greater 
than a certain value, then numerous skips of sizable intensity occur 
during the entry trajectory. Information on the heating rates, total heat 
absorbed, and horizontal range during entry, has been obtained from these 
Z functions and is discussed later. The Z functions in figure 6 . 
could be applied to any planet by noting that cpi for earth is equivalent 
to a value ( pr ) e -1 ^ 2 times as great on another planet, and that a given 
(L/D) for earth is equivalent to a value ( pr ) e -1 / 2 times as great. 

Atmosphere braking for various L/D (ui >1, cp^ < 0).- In entering 

the atmosphere of a planet from space, the approach velocity can be com- 
parable to escape velocity (ui = J2) . It is uneconomical in weight to 
use chemical rockets for reducing the approach velocity in outer space, 
and it is possibly uneconomical in time to use a low- thrust space engine. 
Hence there is considerable interest in the braking process of making 
successive passes through an atmosphere in order to reduce stepwise the 
velocity and the eccentricity of an orbit to near circular conditions 
(ui s 1 ) . 
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In analyzing the atmosphere braking process, cos cp can safely be 
replaced by unity, so the basic differential equation (21) becomes once 
again the same as equation (44) , but the initial conditions are now 

Z(ui) = 0 Z‘(ui) = v/pr 9i 

where uq > 1. By arbitrarily selecting various values for the angle of 
entry cpi , various solutions are obtained corresponding to single passes 
through the atmosphere at various altitudes from the surface. It might 
be more convenient in describing a single pass to select as the arbitrary 
parameter the velocity u ex at the exit of the pass, or 30(uZ)*nax> which, 
for the earth, would be the maximum deceleration in g's experienced during 
this pass, and would be independent of W/CpA. 

In figure 7, four Z functions are presented for nonlifting vehicles 
which start the braking process with essentially escape velocity (ui=1.4) 
but with different values of maximum deceleration in the first pass. The 
short-dash curve (a) corresponds to a maximum deceleration in the first 
pass of 30(uZ) max = 0.46. It is seen that, starting with this initial 
pass (and with no further control exercised on the vehicle) six passes 
would occur before the seventh pass completed the entry process . The 
long-dash curve (b) in figure 7 corresponds to 30(uZ) max = I .65 for the 
first pass. In this case only two passes occur before the third pass 
completes the entry. The other two curves (c) and (d) in figure 7 corre- 
spond to conditions wherein the first pass is the only one, inasmuch as 
it is made sufficiently close to the planet surface to complete entry 
without ever emerging from the atmosphere . 

In computing the Z function for a successive pass , the initial 
angle was assumed to be the same as the exit angle of the previous pass. 
The exit angle was taken at the point where dr/r = du/u. Further dis- 
cussion of these Z functions, and the results of other such functions 
computed for atmosphere braking are presented later. 


RESULTS AND DISCUSSION 


From the various Z functions presented, it is relatively easy to 
study the influence on entry motion of several variables of practical 
interest. For example, we could study the effect of lift-drag ratio on 
deceleration and aerodynamic heating, or the effect of a small error in 
initial angle of descent on the range over which the re-entry process 
takes place. Before considering such topics, however, it is desirable 
to discuss two preliminary items. First, we compare some results from 
the present approximate analysis for an exponential atmosphere with more 
exact machine calculations for a standard atmosphere . This serves to 
provide a feeling for the accuracy of the present analysis, and also to 
show how any of the subsequent results readily can be corrected, if 
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desired, for atmospheric temperature variations. Second, we discuss the 
relative deceleration and aerodynamic heating of various planetary atmos- 
pheres . This provides multiplication factors which enable any of the 
subsequent results for the earth’s atmosphere to be quickly converted to 
results for other planetary atmospheres . 


Comparison of Present Analysis With Other Calculations 


An insight into the approximate accuracy to be expected from the 
present analysis can be obtained by comparison with machine calculations 
of the pair of motion equations for specific vehicles. Differences 
between the present analysis and more exact calculations can arise 
inasmuch as the present analysis makes certain assumptions about the 
trajectory (that is, |dr/r |«|du/u| and |(L/D)tan cp| « 1) and about the 
atmosphere (p^ ~ e~0y) which need not be made in numerical machine cal- 
culations. The a posteriori check of the trajectory assumptions, as 
presented in appendix A, shows that insofar as convective heating and 
peak decelerations are concerned, only a few percent difference should 
be expected for vehicles entering from a satellite orbit. A check of the 
assumption of an exponential atmosphere can be obtained by comparing with 
numerical calculations for some standard atmosphere. In figure 8 a com- 
parison is made of the present analysis with numerical calculations from 
the pair of motion equations using the ARDC model atmosphere. These 
numerical calculations were made by M. W. Rubesin and G. Goodwin using 
equations equivalent to (6) and (7) without discarding any terms. The 
curves in figure 8(a) show close agreement of both the altitude and the 
descent angle as a function of velocity. The curves in figure 8(b) show 
similar agreement of the circumferential distance traveled (As /r ) , and of 
the maximum deceleration (within 6 percent) . This small difference in 
maximum deceleration is believed due primarily to the departure of the 
ARDC atmosphere in certain altitude regions from the idealized exponential 
atmosphere of constant p (l/23,500 ft -1 ). 

As noted earlier, the present analytical method can readily be 
applied using semilocal values of if it is desired to make correc- 

tions to the results in order that they more closely represent some 
standard atmosphere. Corrections also can be made to allow for atmos- 
pheric seasonal variations, or for variations with the earth’s latitude. 

In this sense, analytical results for an exponential atmosphere are actu- 
ally more general for global application than numerical results for any 
single standard atmosphere. This can be seen from the results which 
follow. Let us consider the maximum deceleration_for entry from a 
decaying orbit. This occurs at a velocity near u = 0.43 at which point 
Z = 0.64 (fig. 4(a)). The approximate altitude at which maximum deceler- 
ation occurs is obtained by substituting either into equation (31a) to 
yield (p /p ) , = 5.5(W/C D A)xlO , or into equation ( 31 b) to yield 

0 max du/dt 
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■) 


ft 


(49) 


and is seen to depend on w/CpA. Since N /pr for the ARDC model atmos- 
phere depends weakly on altitude, as shown in figure 2, the maxim-urn decel- 
eration \l 3r(uZ) max = O .278 in a standard atmosphere also wild, 

depend weakly on altitude , and hence weakly on W/CpA. The resulting 
values of -a^/g are shown by the solid curve in figure 9, for (w/CpA) 
values ranging from 0.01 to 1000 lb /ft. They agree very well with the 
points shown which represent numerical integrations (Rubesin-Goodwin) of 
the complete equations for the ARDC model atmosphere. 

If desired, similar corrections for atmospheric variations also 
could be made to other quantities computed for a mean value ,/0r = 30. 
Thus, the distance traveled varies as (|3r)" 1/2 (eq. (28)), and the con- 
vective heating rate varies as (pr) 1,/4 (eq.(36)). It is noted that the 
fluctuations in (3r with altitude, as plotted in figure 2 correspond 
very closely to the fluctuations in T -1 / 2 , as should be expected, since 
|3 = Mg/RT. Hence any variations in mean atmospheric temperature, such as 
seasonal variations or longitudinal variations, can just as readily be 
corrected for as variations with altitude. 

Gazley (ref. 8) has developed an approximate theory for the case of 
orbital decay with L/D = 0 by assuming ucp is constant. This arbitrary 
restriction yields results for orbital decay without lift that are quali- 
tatively similar to the present analysis , but quantitatively dissimilar. 
For example, the density-velocity relationship near peak heating (u = 0.8) 
differs by a factor of roughly 2. For higher velocities the discrepancy 
rapidly increases, and for lower velocities it decreases, me peak decel- 
eration occurs at lower velocities and is not greatly affected by Gazley' s 
assumption. For the earth (*]Jr = 30) he obtains a maximum of 9.6 g at 
u = 0.5k, whereas the present analysis, which does not make any assumption 
about the u(cp) relationship, yields 8.3 g at u = 0.43. 


For a given size and shape of vehicle the deceleration, laminar 
heating rate, total heat absorbed, and Reynolds number vary, according 
to equations ( 26 ), ( 36 ), (39), and (33), respectively, as 


Relative Deceleration, Heating, and Reynolds Numbers 
For Entry Into Various Planetary Atmospheres 


a ~ g Z 

q - PT“ 2/ 3 Mo l/2g3/2 r 5/4plA Z l/2 

Q ~ Pr" 2 / 3 u l/2 gr 5/4 p- 1/4 Z- 1/2 
o 

Re ~ g 1 / 2 p l/2 (i 0 _;L Z 


(50) 
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In the case of nonlifting entry from a decaying orbit (zero initial angle 
of descent) , the characteristics of the planetary atmosphere CTO do not 
enter the differential equation or the initial conditions (Zi = 0, 

Zi’ = 0); hence Z can be disregarded in computing the relative values 
of the above quantities for various planets . In the mor e ge neral case of 
entry from high altitude (Zi = 0) with fixed values of 3r cp j_ and pr ( L /D ) , 
the Z function still would be the same for all planets. Neglecting dif- 
ferences in Prandtl number and ratio of specific heats, we have for several 
planets the following relative values applicable to nonlifting entry fro m 
decaying orbits, or to any other type of entry where the values of cpj_ 

and -v/p r(L/D) are fixed: 



(du/dt) e 

relative 

deceler- 

ation, 

(g -s/Pr)® 

relative 
heating rate, 

(P 0 l/2 g 3/2 r 5/4 P 1/ 4 )® 

relative total 
heat absorbed, 

(p 0 1 ^ 2 gr 5/4 P _1/ 4 ) e 

relative 
Reynolds 
number , 

Venus 

0-9 

0.7 

0.8 

1 . 

Earth 

1.00 

1.00 

1.00 

1.00 

Mars 

.2 

.09 

.2 

.4 

Jupiter 

5. 

50. 

50. 

2. 


It is to be remembered that, in the case of a vehicle with lift, in order 
to have the same Z function, a given L/D ratio on earth is equivalent 
to an L/D ratio (pr) e " l/2 times as great on a planet other than earth, 
and that a given (pi on earth also is equivalent to a cpi value 
(3r)®" l/2 times as great on another planet. This equivalence, together 
with the above table, enables any result for earth to be converted to a 
result for each of the other planets. 

In the special case of entry at a constant angle 9i, the atmospheric 
characteristics enter the initial conditions on Z (Zi* = s/f3r sin cpi) . 
Since equation (4o) shows that Z ~ ,/jjr for this type of re-entry, we 
include this in the expressions (50) to obtain the following relative 
values applicable only to ballistic entry (L/D = 0) at constant cp: 



(du/dt)® 

relative 

deceler- 

ation, 

(gpr) e 

relative 
heating rate , 
(Po l/2 g 3/2 r 3/2 P l/2 )e 

Qe 

relative 
total heat 
absorbed, 
(po l/2 S r )® 

Re® 

relative 

Reynolds number, 
(g 1/2 3r 1/2 ia 0 " 1 ) 3> 

Venus 

0.9 

0.7 

0.7 

1 . 

Earth 

1.00 

1.00 

1.00 

1.00 

Mars 

.09 

.06 

.2 

.2 

Jupiter 

n. 

70. 

20. 

4. 
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These relative values for ballistic entry are exactly the same, of course, 
as would be obtained directly from the theory of reference 1 and are 
applicable for initial angles greater than about 5°. The previous table 
would apply for = 0°. For nonlifting entry with cpi the order of a few 
degrees, the relative values for various planets would be intermediate to 
the above two tables. 

We see that entry into the atmosphere of Venus involves only slightly 
less deceleration and heating than does entry into the earth’s atmosphere, 
whereas entry into Mars involves much less deceleration and heating, and 
entry into Jupiter, much more. The Reynolds numbers, however, are not as 
greatly different for the various planets. 


Effect of Lift on Deceleration, Heating Rate, and Total Heat 
Absorbed During Entry From Decaying Orbits 


From the relative values of deceleration and heating for various 
planets, together with the Z functions already presented, certain 
quantities of practical interest readily can be computed. The remainder 
of this report concerns such application of the Z functions for the 
various types of entry. In the present section we discuss first lifting 
entry from decaying orbits (uq = 1 , cpi = 0 ) . 


Deceleration.- A plot of the horizontal deceleration du/dt in g’s 
for the earth's atmosphere (equal to 30 uZ) is presented in figure 10 as 
a function of the dimensionless velocity u for various lift-drag ratios. 
The powerful effect of L/D ratios the order of only a few tenths is 
evident from this figure. It is also evident that the maximum decelera- 
tion occurs near a velocity of u = 0.4. These curves are independent of 
the shape, size , and mass of the vehicle. The re sultant deceleration is 

taken as a = J (du/dt ) 2 + [(dv/dt) - (u 2 /r) + g] 2 . For no motion this 
expression reduces to g, the gravitational constant of the planet. By 
substituting equations (l 6 ) and ( 20 ) we have 


a 

g 


n/Pt uZ 
cos cp 




(51) 


which, for small angles ( |cp| « L/D, cos cp £ 1, tan 2 cp « 1) yields 



•ypr(HZ) max ^ + ( L / D ) 2 


A plot of this approximation for the maximum resultant deceleration is 
shown in figure 11 for the several planets considered. Once again the 
strong influence of the L/D ratio near L/D = 0 is evident. Also 
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evident are the relatively low decelerations for Mars compared to earth 
and Venus, and the relatively high values for Jupiter. 

From the viewpoint of human tolerance to acceleration stress, it is 
not only the peak deceleration which must be considered, but also the 
orientation of the body, the duration of stress, and the rate of onset of 
deceleration. Numerous experiments with the human centrifuge have shown 
that h uman tolerance is greatest in transverse orientation; that is, with 
either chest-to-back or back-to-chest loading. Centrifuge experiments 
(see, e.g., ref. 16 and the references quoted therein) also have shown 
that the magnitude of acceleration is relatively more important than the 
duration, in the sense that if the acceleration is increased 10 percent, 
the tolerable duration is decreased by a factor of about 2. Thus, a 
method believed to be conservative for calculating the effective duration 
At during entry is to assume that the maximum deceleration acts over the 
entire time it would take for this deceleration to slow the vehicle from 
orbital velocity to rest. Curves of maximum deceleration versus duration 
computed in this manner are presented in figure 12 for various planetary 
atmospheres and for various L/D ratios. Included in this figure is a 
boundary representing human tolerance in the transverse orientation for 
conditions of rapid onset of acceleration (ref. 1 6 and references quoted 
therein) . This boundary also is conservative inasmuch as entry decelera- 
tions are built up relatively slowly under which conditions, according to 
the centrifuge experiments of reference 17 , the body circulation builds 
up a reflex action of effectiveness comparable to that provided by a 
G-suit. The conservative limits determined from this figure are indicated 
in figure 11. It is evident from both figures 11 and 12 that the deceler- 
ations for orbital entry into the earth’s atmosphere are well within human 
tolerance even for nonlifting bodies. For Mars, human tolerance is suf- 
ficient to permit entry at sizable angles of descent or with negative 
lifting devices. Manned entry into Jupiter, however, would require a 
positive lifting body, or some other device in order to maintain the 
decelerations within human tolerance . 

Heating rate.- In examining the effect of lift on convection aero- 
dynamic heating of entering vehicles, we can use the same Z functions 
as employed in studying the decelerations. We note first that for many 
vehicles, the values of Reynolds number near peak heating are sufficiently 
low that one would expect a considerable extent of laminar flow, yet suf- 
ficiently high to be in continuum-gas flow rather than f ree-molecule flow. 

A plot of Re/Z at peak heating is presented in figure 13 as a function 
of w/CpA for entry from orbital decay into the earth’s atmosphere. A 
vehicle on a large parachute would correspond to W/CpA the order of 
0.1 lb /ft 2 , and, with L/D = 0, to Re of about 10 3 ft -1 . For such con- 
ditions the peak heating, which occurs at a Mach number Moo = 20, would 
be near the slip-flow regime (Re / m^ = l) . A reasonably blunt metallic 
structure would correspond to W/CpA values the order of 10 to 100 lb /ft 2 , 
and to values of Re/z the order of 10 3 to 10 5 . Such values are well 
within the continuum regime, yet low enough to be associated with laminar 
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flow. The curves in figure 13 are for earth hut can he applied to other 
planets hy multiplying the ordinate hy the value of the relative Reynolds 
number already tabulated for several planets . 

For a given atmosphere the laminar heating rate is proportional to 



A plot of the dimensionless heating rate q as a function of u is 
presented in figure lb for entry from decaying orbits. The maximum value 
occurs at a velocity u of about 0.8 and is a function only of the param- 
eter n/ ( 3r (L/D) as follows: 


•/Ml/d) 

( L//D ) Earth 

3-max 

77 1/4 

3-max 

-15 

-0.5 

0.375 

0.783 

- 7-5 

-.25 

.302 

.741 

-3 

-.1 

.253 

.709 

0 

0 

.218 

.683 

3 

.1 

.184 

.656 

7-5 

.25 

.138 

.610 

15 

.5 

.098 

.560 

30 

1.0 

.070 

.514 


For l/d ratios greater than 1 the asymptotic solution Zjj = — - k_: _i 

N /P r u(L/D) 


can be used to yield q max = 


J 3 *s/3r(L/D) 


as noted in appendix C. 


We will consider that the vehicle dimensions and weight (R, A, and W) 
are fixed, and will study the influence of vehicle shape (Cp and l/d) . 
Under these conditions the maximum heating rate is proportional to 

. The effect of lift-drag ratio on maximum heating rate (which 
occurs at a u of roughly 0.8) is illustrated in figure 15 for entry 
from decaying orbits. The quantities plotted have been normalized to 
unity for L/D = 0, and can be applied directly to any planet, as can the 
curves in figure l4, by recalling that a given L/D for Earth is equiva- 
lent to a value ( 3r ) times as much for a planet other than Earth. If 
the L/D ratio could be increased indefinitely without changing the drag 
coefficient such as by using reaction lift, then the maximum laminar heat- 
ing rate would be proportional to the dotted line in figure 15 represent- 
ing and would decrease indefinitely with an increase in L/D 
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(asymptotically as (l/d)" 1/2 for L/D greater than about 0.5) • Physically 
this decrease arises because the greater the lift, the less rapid the vehi- 
cle descends, so that the heating occurs at higher altitudes where the 
density is lower. On a practical device which uses aerodynamic lift, how- 
ever, the L/D ratio cannot be increased much without making the vehicle 
more slender and decreasing Cp; a decrease in Cp increases the heating 
rate (~ 1 /Jcj)) because it results in less slowing down, thereby causing 
the peak heating to occur at lower altitudes where the density is higher. 

As a result, there is an optimum l/d ratio for minimizing the heating 
rate which, for the three families of shapes indicated in figure 15, is 
near the range of L/D between about 0.5 and 1. For the family of half- 
cones and half -paraboloids, the l/d ratio was changed by changing the 
fineness ratio while maintaining the flat tops parallel to the stream 
direction. For the family of flat plates the l/d ratio was changed by 
changing the angle of attack. In all cases, Cp and l/d were computed for 
Newtonian flow. The optimum L/D ratio is seen to depend somewhat on the 
particular aerodynamic shape, since L/D and Cp are coupled somewhat dif- 
ferently for different shapes. It is evident that the net benefit to be 
gained by using aerodynamic lift amounts to about a factor of 2 in reduc- 
ing the maximum rate of aerodynamic heating at a stagnation point. 

Inasmuch as the optimum l/d ratios for minimizing the maximum heat- 
ing rate are greater than about 0.5, they are in the range where the Zji 
function for orbital decay is a good approximation near peak heating (see 
fig. 4(b)). From equations ( 36 ) and (4l) we see that for a given planet 
and given radius at a stagnation point. 


and we see that the various minima in figure 15 each correspond to enter- 
ing at C]^ aY . The peak heating always occurs at a dimensionless velocity 

u = J 2 J 3 = 0.82. For flat plates in Newtonian flow C] ^ 3Y = 0.77 at an a 

angle of attack of 55°, for which L/D = 0.71. As noted in appendix^, 
these conditions also turn out to represent optimum ones for minimizing 
the total heat absorbed for skipping-type entry, because in this case also 
q and Q vary as (l/d)" 1/2 . 

Surface temperature for radiation equilibrium .- The stagnation sur- 
face temperature experienced during entry of a structure having relatively 
small heat capacity (e.g., a thin skin) is calculated by equating the 



since L/D = Cl/Cp, 



(54) 
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radiation heating rate to the convective heating rate . For entry from 
decaying orbits we may set cos cp = 1, inasmuch as cp near peak heating 
varies from -2.6° to -0.2° as L/D varies from 0 to 1. We have 

ecrT Ws 4 = q s = 590 1 Btu ft“ 2 sec _1 (55) 

where e is the surface radiative emissivity, R the radius of curvature, 
and 0 = 0.48xl0 -12 Btu ft -2 sec -1 °R -4 is the Stefan-Boltzmann constant. 

By substituting the value of g c for Earth there results (for T v in °R, 
R in ft, w/CpA. in lb ft" 2 ), 

T„ s € 1,4 R 1,s = 3840 fe) 178 ? 1 ' 4 (56) 

_ 5 / 2 

where q for laminar flow is equal to u Z 1/2 . The maximum value 
q mAy 1/4 is listed in the preceding table for entry from decaying orbits. 

Other types of entry would require the use of other Z functions, but 
equation (56) would remain unchanged. For a planet other than Earth, the 
radiation-equilibrium temperature calculated from the above equation for 
Earth would have to be multiplied by the l/4-root of a quantity already 
tabulated; namely, the relative rate of heating q e for that planet. The 
relative radiation-equilibrium temperature factors T w ^ = q$ 1/4 are: 



Venus 0 . 91 
Earth 1 . 00 
Mars . 55 
Jupiter 2 . 7 


A graph of the maximum temperature parameter Tv s e 1/4 R 1/8 for entry from 
decaying orbits is presented in figure 16 as a function of W/CpA (W in 
Earth weight) . It is noted that the numerical calculations for nonlifting 
satellites descending in the Earth's atmosphere, as reported by Kemp and 
Riddell (ref. 6) and by Gazley and Masson (ref. 5 ) > agree well with the 
analytical variation represented by the present analysis. 

The curves for T Wg in figure 1 6 could be applied to other planets 
for any given value of ./pr (L/D) by multiplying the ordinate by the 
quantity T w ^ tabulated above. Since L/D is a more convenient variable 
than v/pr (L/D), however, a separate plot of the parameter 

T Ws e 1/4 /(W/CpA R) 1/8 (which represents the maximum surface temperature that 
is experienced during entry for radiation equilibrium at a stagnation 
point of radius of curvature R) is presented in figure 17 as a function 
of L/D for several planets. The coupling between Cp and l/d is taken 
as that for the family of half -paraboloids . The others would not be 
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greatly different, as may be seen from the curves in figure 15* We can 
deduce from figure 17 , for example, that a nonlifting body with e = 0.9; 
and w/CdAR = 1 lb ft' 3 (e.g., R = 1 ft and W/CpA = 1 lb ft" 2 or R = 10 ft 
and W/CdA = 10 lb ft' 2 ) would experience during entry from orbital decay 
a maximum stagnation temperature of approximately 1000° F for Mars, 

2000° F for Venus, 2200° F for Earth, and 6800° F for Jupiter. 

Total heat absorbed.- It is emphasized that the effect of lift-drag 
ratio on the total heat absorbed Q is quite different from the effect 
just discussed on the heating rate q. The use of lift prevents a vehicle 
with a given drag coefficient from descending as rapidly as a nonlifting 
one, thus leading to lower heating rates at higher altitudes, but the lift 
also prolongs the descent markedly. This prolongation dominates over the 
reduced rate of heating, to lead to a net increase in total heat absorbed 
with increasing L/d. That the total heat absorbed must increase with an 
increase in L/D, may be clearly seen from the general equation 


developed by Allen and Eggers in reference 1. For a given Cp, an increase 
in l/d does not change the kinetic energy loss, but it does increase the 
effective laminar skin-friction coefficient Cp' inasmuch as the corre- 
sponding increase in altitude results in the heat being taken aboard at 
lower Reynolds numbers where Cp ’ is higher. 

The quantitative magnitude of the increase in Q, with an increase 
in L/D may be deduced from equation (39) for Q (which neglects the heat 
radiated from the surface). For a given atmosphere (given Pr, p, g, r, 
p) and a given size and weight (A, R, W), Q for laminar flow and 
cos cp = 1 is proportional to the quantity 


where Q is a function of l/d and is very insensitive to the_lower 
limit u down to which the integration is carried (providing u is 
small) . For convenience in evaluating Q from the Z functions, we 
select an arbitrary upper limit ui = 0.99* The following values for 
Q are obtained for entry from decaying orbits: 



(57) 
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*7pF(l/d) 

( L / D )Earth 

Q 

for ui = 0.99 

-30 

-1.0 

0.75 

-15 

-.5 

•93 

-7-5 

-.25 

1.09 

-3 

-.1 

1.23 

0 

0 

1.36 

3 

.1 

1.5^ 

7.5 

• 25 

1.90 

15 

.5 

2.53 

30 

1.0 

3-5^ 


For L/D ratios greater th an 1 the asymptotic Zjj function can he used 
to yield Qjj = O.62(0r) 1/4 ^/l/d for the heat absorbed between u x = 0*99 
and u = 0. (See appendix C for a more general expression for Qji . ) 

The effect of lift-drag ratio on the total convective heat absorbed 
(disregarding radiation from the surface) during entry into the earth’s 
atmosphere from decaying orbits is plotted in figure 18. These curves 
are normalized to unity for l/d =0. In contradistinction to the_effect 
on q, an increase in l/d by itself is seen to always increase Q, and 
hence Q, as anticipated from equation (57) • When the coupling between 
l/d and Cp is considered, an optimum occurs at negative l/d ratios, 
near the range -0.7 to -0.5 . In view of the fact that these negative 
l/d ratios result in high decelerations (fig. ll) they would not be 
feasible for a manned entry into the earth’s atmosphere; the practical 
optimum for a heat-sink vehicle would be near L/d = 0. 

In figure 19 curves are presented of the total heat absorbed per unit 
area during entry into various planets from decaying orbits . Radiation 
from the surface is disregarded for these curves . They represent the 
family of half -paraboloids , but the other families would not be signifi- 
cantly different. As would be expected, the minimum for each planet occurs 
at a negative L/D ratio. For Mars the decelerations are not excessive 
for L/D near -0.5 (see fig. 11) but the reduction in total heat absorbed 
compared to a nonlifting vehicle is only about 10 percent . 


Nonlifting Entry From Deflected Orbits 


In the discussions thus far we have considered only the trajectories 
resulting from decaying orbits wherein the initial descent angle is 
essentially zero. This type of entry leads to relatively shallow angles 
of descent with relatively low heating rates, but provides very little 
control over the time of entry and the location of impact . One method 
commonly envisioned to fix the time of entry, and greatly improve the 
accuracy of landing in a predetermined area, is to induce entry by 
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suddenly deflecting an orbit so as to enter at some initial flight path 
angle cpq . A retrorocket force , for example , or a rocket force applied 
in the direction toward the planet center will initiate such entry. 

Induced entry of this type, however, results in greater decelerations and 
can affect the aerodynamic heating problem either adversely or favorably. 

A curve is presented in figure 20(a) showing the effect of initial angle 
cpj_ on the maximum deceleration experienced during entry of nonlifting 
vehicles into the earth's atmosphere. Also shown for comparison is the 
approximate limit of human tolerance (for rapid onset with transverse 
orientation) , and a dotted curve corresponding to the Allen-Eggers theory 
for cp = constant = cpp. This theory for ui = 1 can be used for descent 
angles greater than about 4° or 5° • Above about -cpj_ = 3° "the decelera- 
tions exceed human tolerance, so that some method of deceleration allevia- 
tion, such as provided by lift, or by increasing the value of W/CpA 
during descent, would have to be employed for manned vehicles entering at 
these larger angles of descent. The curve of (du/ dt)mav in Earth g*s 
can be applied to any planet by regarding the abscissa scale as being 
-(s/p r)_cpi and then multiplying the ordinate scale by (g pr) e . 

The effect of initial angle on maximum laminar heating rate and on 
the total laminar heat absorbed is shown in figure 20(b) . As would be 
expected, the steeper the descent the greater the heating rate. The total 
heat absorbed, however, is less for the steeper descents because the 
shorter duration more than compensates for the greater laminar heating 
rates. Equation (57) shows that this must be the case, since entry at 
larger angles results in the heat being taken aboard at lower altitudes 
where the laminar skin-friction coefficients are small. If the flow were 

turbulent the corresponding reduction in Cp * and hence in Q with an 

increase in descent angle would be less. The curves in figure 20(b) 
approach the curves developed from the Zj function corresponding to the 
solution of Allen and Eggers (see eqs. (C3) and (c4) of appendix C). In 
order to be consistent with the other values of Q representing the heat 
absorbed from u = 0.99 "to u = 0, a calculated factor 0.8A has been applied 
to equation (cU) which represents the heat absorbed from u = 1 to u = 0. 

It is seen from figure 20(b) that the Allen-Eggers solution for heat 

transfer in this case (uj[ = l) is quite accurate for descent angles greater 
than about 2°. The curves in figure 20(b) can be applied to other planets 
by regarding the abscissa as a scale for the quantity -(sffir) - 

In the figure 20 (c) a curve is presented showing the strong influence 
of initial descent angle on entry range for Earth. Two incremental ranges 
are shown: a solid line curve for the distance between the point where 

u = 0.995 and the impact point (u = 0) , and a dashed-line curve for the 
distance between u = 0.99 and impact. From the slope of the solid-line 
curve we obtain the lower curve shown of average miss distance for an 
error in cpj_ of 0.5°* It is to be remembered that this miss distance 
curve does not consider the essentially dragless portion of a deflected 
orbit from the point of orbit deflection to the point where u = 0.995* 
and hence it is indicative of only the entry portion of the practical 
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problem of estimating miss distance. The curve illustrates, however, the 
advantage of using a small initial descent angle in order to greatly 
improve the ability to determine impact point. 

A further contribution to miss distance which can be studied with 
the present equations is that due to atmospheric variations in temperature 
with either season or latitude. Equation (28) shows that As 1 , 

so that a ±15-percent seasonal variation in temperature would correspond 
to a +7-percent variation in ,/pr and in As. For small initial angles, 
say cpi = -1°, the range during entry from u = 0.995 to impact is roughly 
1000 miles according to figure 20 (c), and hence the impact point would 
vary ±70 miles. The entry range would be greater in summer than in winter. 

A graph of the Reynolds number per foot at peak heating for nonlifting 
entry into the earth’s atmosphere with Vi = 1 is presented in figure 21 
for -cpi = 0°, 5°, 10°, 20°, 40°, and 90°. The -9i = 0° curve is based 
on the Z function of f igure 4(a) . All others are based on the Zj 
function corresponding to the Allen-Eggers solution. Entry at other values 
of Vi, according to this solution, results in values of Re proportional 
to Vi. 


Lifting Entry From Deflected Orbits 


If a vehicle with l/d > 0 enters the atmosphere from a deflected 
orbit at a sufficiently large initial angle of descent, the entry trajec- 
tory is comprised of one or more skips. This is to be expected on physical 
grounds and is evident from the Z functions already presented in fig- 
ures 6(b) to 6(e). During the first portion of descent, a vehicle under- 
going a sizable skip will, at the bottom of the skip, decelerate and take 
on heat at a lower altitude than a vehicle at the same velocity which 
glides in smoothly from a decaying orbit (cpi = 0) • For large initial 
angles of descent, then, we might expect a skipping vehicle entering from 
a deflected orbit to experience greater decelerations, higher heating 
rates, and shorter entry range than a gliding vehicle entering from a 
decaying orbit. On the other hand, since the skipping vehicle takes on 
most of its heat at a lower altitude (where the skin-friction coefficients 
are lower) we would expect from equation (57) that the skipping vehicle 
would absorb less total heat during entry than the orbiting-decay vehicle. 
Calculations from the Z functions of figures 6(b) to 6(e) show these 
various expectations to be the case for initial descent angles -(TiJEarth 

greater than about 1°. This is illustrated in figure 22(a) for maximum 
laminar heating rate, in figure 22(b) for total laminar heat absorbed, and 
in figure 22(c) for entry range. The expected increase in deceleration is 
already evident from figures 6(b) to 6(e) which show 30 uZ ~ du/dt as the 
ordinate . 

If a vehicle with L/D > 0 enters the atmosphere from a deflected 
orbit at a very small initial angle of descent, so that the trajectory 
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might be described more appropriately as a rippling descent rather than 
a skipping one, then the peak deceleration and maximum heating rates can 
actually be slightly smaller than for the same vehicle gliding in from 
a decaying orbit. What happens in such cases may be seen, for example, 
in figure 6(c) by comparing the curves for -q)j_ = 1° and -cpj_^= 0°. The 
rippling entry (-cpi = 1°) has one maximum on each side of the maximum 
for -cpj_ = 0° representing orbital-decay entry. These two maxima in 
deceleration for — cp-j_ = 1° are slightly less than the single maximum for 
-cpj_ = 0°. A similar situation can exist for the maxima in heating rate. 
As a result, the curves in figure 22(a) for the dimensionless maximum 
heating rate q max for lifting vehicles entering from deflected orbits 
show slight waviness and sometimes slight reductions below the values for 
-cpj[ = 0 when the initial descent angle is less than about l/2 to 1°. 
Consequently, we can say that, in principle, a rippling-type descent from 
a deflected orbit can have lower maximum heating rates than a gliding 
descent, but for practical purposes, there is no significant difference 
between the two. 


Composite Entry 


It may be desirable to combine lifting and nonlifting entry in order 
to achieve some advantages of both types. For landing maneuverability it 
obviously is advantageous to employ a lifting vehicle. The total heat 
absorbed by a lifting vehicle, however, is much higher than for a nonlift- 
ing vehicle (fig. 18) . The optimum use of aerodynamic lift reduces the 
maximum heating rate only to about one-half that of a nonlifting vehicle 
of the same W/A. Nonlifting vehicles can more easily be constructed 
with much lighter W/A ratios by employing, for example, a large, light 
drag device (for example, a parachute). The larger the device, the smaller 
is the heating rate (q ~ l/ >JAI ~ , the smaller the entry Reynolds 

numbers (Re ~ (w/CpA )l ~ 7" 1 ), and the better the possibilities are of 
maintaining laminar flow. Nonlifting vehicles with shuttlecock stability 
are advantageous also from the viewpoint of minimum control requirements 
during entry. Hence, an evident composite type of entry, which combines 
some of the desirable features of lifting and nonlifting trajectories, 
would be to enter first without lift but with a small w/CpA provided 
by a drag device; then, when the velocity is reduced to a certain value 
u-b the device is jettisoned or retracted, leaving a lifting vehicle of 
larger W/CpA for the remainder of the descent. 

A practical compromise is required in selecting u^, because the 
drag device should be jettisoned as soon as possible from the viewpoint 
of achieving maximum maneuvering range, but as late as possible from the 
viewpoint of achieving major reductions in heating rate. For the initial 
nonlifting portion of descent let the drag-weight parameter be (w/CpA) 0 
and the Z function be Z Q . For the subsequent portion let the corre- 
sponding quantities be (w/CpA ) 1 and Zi. Since the altitude y and the 
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angle of descent cp are continuous at the break velocity ub, we have 
two conditions from equations (l4) and ( 17 ) 



for determining the initial conditions Zi^ = Zi b and Zi ^ * = Zi^ r for 

the second portion of descent. Hence the Zi function can be determined 
approximately from equation (C13) of appendix C by substituting up = u^, 

Z± = sin 9i = 9b> and cos 9 = 1 . The maximum heating rate occurs 

near the bottom of the first dip after the break, and can be obtained 
from equation (C24) in appendix C with the same substitutions. The total 
heat absorbed in this dip can be obtained from equation (C25) . 

As an example let us consider the case of a large drag device 
((w/CdA)i » (w/CdA) 0 ) jettisoned at a velocity ub during entry from 
a decaying orbit. In order to minimize the peak heating after jettisoning, 
as well as minimize the total heat absorbed during the skip, a value 
l/d = 0.7 is selected. Curves showing the resulting values for maximum 
heating rate q^ after jettisoning, and total heat absorbed Qx during 
the first skip, are presented in figure 23 as a function of the break 
velocity u^. We see that a large drag device carried down to U 5 = 0.4, 
for example, would have a maximum heating rate about l/4 of that for the 
same vehicle with no drag device. 

It may be noted that the deceleration history for a drag device 
jettisoned at 0^ = 0.4, for example, is essentially the same as the 
acceleration history investigated in the human centrifuge tests of refer- 
ence 18. The select individuals for these centrifuge tests did not 
blackout (or grayout, or even get dizzy) during the runs. They were able 
to perform continually simple dual control operations even when the accel- 
eration dropped suddenly from about 8 g to about 2 g. 


Comparison of Several Types of Entry With u^ = 1 


It is interesting to compare the relative magnitude of aerodynamic 
heating for the several types of entry discussed. The dimensionless 
maximum heating rate q^ and the dimensionless total heat absorbed Q 
are used for this comparison. They would be proportional to the actual 
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heating rate and the total heat absorbed for vehicles of the same size 
and w/Cj)A. The table which follows summarizes these quantities for 
seven different types of entry, all starting with Uj_ = 1. 


Type of entry 

L/D 

^max for 

O 

II 

Near optimum glide, for 
minimum q^y = 0) 

0.7 

0.084 

3.0 

Near optimum ripple for 

minimum q maY (-cp i = 0.5°) 

• 7 

.083 

2.9 

Near optimum glide, for 
minimum Q(cpj_ = 0) 

- .5 

.78 

.93 

Near optimum first skip 
for minimum Q 

• 7 

•15(-<Pi = 2°) 

.90 

Nonlifting (cp ± = 0) 

0 

.22 

1.4 

Nonlifting, from deflected 
orbit with -cpi = 2° 

0 

.27 

.93 

Composite, large drag device 
jettisoned at u-^ = 0.4 

0 for u > 0.4 
0.7 for u < 0.4 

.02 

.16 


In comparing these values it should be remembered that_the actual quanti- 
ties of interest for a given w/A are Qmax/ Q/ an< ^- that 

nonlifting vehicles are placed at a small disadvantage in the table 
because they presumably can be designed with somewhat higher values of 
C D than lifting vehicles. It is noted that the total heat absorbed in 
the case of the skip vehicle, corresponds only to the first skip. Pre- 
sumably this is all that should be considered if the vehicle is designed, 
as suggested by Ferri (ref. 7), to radiate essentially all of the heat 
absorbed after each skip. 


Atmosphere Braking 


During entry of a planet's atmosphere from space at near escape 
velocity, possibly severe deceleration and heating problems can occur 
during the process of passing through an outer segment of the atmosphere. 
The closer a pass is made to a planet surface, the greater is the braking 
action, the greater the deceleration, and the greater the rate of aero- 
dynamic heating. The Z functions for four different entry histories of 
nonlifting vehicles starting with escape velocity (ui = 1.4) have already 
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been presented in figure 7- These functions apply to any planet. They 
are based on the assumption that after the initial pass no further control 
of the vehicle is exercised. 

Entry (a) is initiated with 30(uZ) max = 0.46 during the first pass 
(0.46 g maximum deceleration for earth) and corresponds to a dimensionless 
peak heating rate of q mRy = 0.24 at u = I. 38 . The successive peaks 
correspond to qinax progressively less, while the seventh pass, which 
starts from u = 1.08 and completes the entry, corresponds to q^g^. = 0 . 20 . 
As might be expected this is not far from the value 0.22 corresponding to 
orbital decay from = 1 with l/d = 0. Since is a measure of 

the maximum temperature experienced by a radiation-cooled vehicle, it 
follows that entry of such a vehicle could be completed on the seventh 
pass, without the temperature during any of the atmosphere braking passes 
exceeding appreciably that experienced during orbital decay. 

Entry (b) in figure 7 is initiated with SOCuZ)^^; = I .65 in the first 
pass during which an amount of heat is absorbed corresponding to Q = 1 . 5 * 
This heat could be radiated to space before the second pass is made in 
which an additional amount Q = 1.4 is absorbed. The third pass starts 
from u = 1.09 and completes the entry with Q = 1 . 7 . These values are 
not far from the value Q = 1.4 corresponding to orbital decay with 
L/D = 0. Since Q is a measure of the total heat absorbed by a heat- 
sink vehicle, it follows that such a vehicle could complete an entry on 
the third pass without absorbing much more heat during each of the two 
atmosphere braking passes than that absorbed during orbital decay. 

Entries (c) and (d) in figure 7 sire completed in a single pass and 
both lose an amount of kinetic energy _(l/2)m(l. 4 ^Tgr ) 2 = mgr. They absorb 
a quantity of heat corresponding to Q = 2.9 and Q = 2 . 1 , respectively, 
and experience maximum heating rates corresponding to q max = 0*58 and 
Imax = 0.73* respectively. The total laminar heat absorbed by (d) is 
less than (c), even though the maximum heating rate is greater, because 
entry (d) corresponds to a closer pass to the planet surface for which 
the heat is taken aboard, on the average, at lower altitudes where the 
friction coefficients are lower (see eq. ( 57 ))* 

In addition to the four Z functions just discussed, a number of 
Z functions (not presented) have been computed for lifting vehicles 
undergoing single atmosphere braking passes in which the entering velocity 
is uq and the exit velocity is u ex * Results are presented in figure 24 
for ui = 1.4 and in figure 25 for uq = 1.2. In each figure curves are 
presented for the maximum value of horizontal deceleration 30(uZ) max , the 
dimensionless maximum laminar heating rate q^^, and the dimensionless 
laminar heat absorbed Q during the single pass. The curves are labeled 
as to the l/d values corresponding to earth; they also can be applied 
to other planets by recalling that a given value of L/D on Earth is 
equivalent to a value (pr) 0 ~ 1 / 2 times as much on another planet. 
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An interesting feature of these results for single atmosphere brakings 
is that for a given loss in kinetic energy (given u e x) > they exhibit the 
opposite variation with L/D from that previously found for orbital 
decay. Thus, an increase in L/D decreases the maximum deceleration for 
orbital decay but increases it for atmosphere braking; an increase in 
L/d decreases the heating rate Q.max for orbital decay but increases it 
for atmosphere braking; an increase in l/d increases the heat absorbed 
Q for orbital decay but decreases it for atmosphere braking. From a 
mathematical viewpoint the reason for this contrasting behavior is that 
the gravity minus centrifugal force term (l - u 2 )/uZ in the basic differ- 
ential equation changes algebraic sign at u = 1. From a physical view- 
point, the effect of L/d on atmosphere braking can be understood by 
noting that in order to lose the same amount of kinetic energy , a lifting 
vehicle must pass closer to the surface than a nonlifting one. Hence at 
the lower altitude the deceleration and rate of heating of the lifting 
vehicle are greater, while the friction coefficients are smaller and hence 
the heat absorbed for a given loss in kinetic energy is smaller (see 

eq. ( 57 )) • 

A plot of the maximum surface temperature parameter 
T WQ e 1 / 4 /(w/CpAR ) 8 as a function of the maximum deceleration in Earth 

g • s is presented in figure 2 6 for atmosphere braking in various planets 
with l/d = 0. These curves are for a single pass starting with uj = 1*4 
It is seen that in the earth* s atmosphere, for example, the maximum decel- 
eration that can be experienced in a single pass and still enable the 
vehicle to exit from the atmosphere at some velocity u e x > 1 , is about 
3.5g. If the nonlifting vehicle attempts to decelerate more than this 
by passing closer to the surface, then before it exits from the atmosphere 
the velocity is reduced to u = 1 at some point within the atmosphere and 
the vehicle completes entry in a single pass experiencing at least 7 * 2 g 
deceleration in the process . Any pass still closer to the surface only 
increases further the maximum deceleration and temperature. When the 
maximum deceleration during a single pass jumps discontinuously from 3-5g 
to 7 . 2 g, the corresponding maximum temperature does not jump because the 
maximum temperature already has been experienced before u = 1 was reached 
The limiting maximum deceleration for atmosphere braking in Mars is seen 
to be much less (0-7 Earth g) , and for Jupiter much more than for Earth. 

A companion plot to figure 26, only for the laminar heat absorbed 
per unit area in a single pass, is presented in figure 27 . These curves 
also are for L/D = 0 and ui = 1.4. In this case, the heat absorbed 
increases discontinuously when the maximum deceleration increases discon- 
tinuously (from 3.5g to 7.2g for Earth) because of the additional loss 
in kinetic energy. Any pass still closer to the surface increases the 
deceleration but decreases the laminar heat absorbed. This decrease 
exists because, for a given loss in kinetic energy, any pass taking on 
its heat at lower altitudes will have smaller laminar friction coeffi- 
cients, and hence less total heat absorbed (see eq.. ( 57 ))* 
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CONCLUDING REMARKS 


An approximate analytical solution for the motion and aerodynamic 
heating of a lifting vehicle entering a planetary atmosphere has been 
obtained by disregarding two relatively small terms in the complete motion 
equations, and then introducing a mathematical transformation which reduces 
the pair of motion equations to a single, ordinary, nonlinear differential 
equation. Relatively few solutions to this differential equation provide 
quite general results inasmuch as the basic equation is independent of 
the physical characteristics of a vehicle , as well as independent of the 
sea-level characteristics of an atmosphere . The solutions apply to any 
exponential planetary atmosphere . 

Certain asymptotic solutions in closed form result from a process of 
truncating various combinations of terms from the basic nonlinear differ- 
ential equation. The aggregate of terms represents vertical acceleration, 
vertical component of drag force, gravity force, centrifugal force, and 
aerodynamic lift force. This truncation procedure yields an asymptotic 
solution for ballistic vehicles entering at relatively steep angles of 
descent (which solution is identical to that of Allen and Eggers), an 
asymptotic solution for glide vehicles of relatively large lift-drag ratio, 
and a solution for skip vehicles. 

Comparison of the present solution for an idealized exponential 
atmosphere with digital computing-machine results for a standard atmos- 
phere reveals differences the order of about ±10 percent. These rela- 
tively small differences are due primarily to the variations in atmospheric 
temperature with altitude in the standard atmosphere. The present analytic 
solution enables corrections readily to be made in order to yield results 
applicable to any standard atmosphere, or to an atmosphere which has 
variations in temperature with season or with latitude . 

Maximum deceleration during entry into an exponential atmosphere 
from a decaying orbit does not depend on the vehicle weight , shape , or 
dimensions; it occurs at a velocity of about OA of orbital velocity, 
and is much less for lift-drag ratios as small as a few tenths than for 
a lift-drag ratio of zero. Even for nonlifting vehicles, though, the 
decelerations are within human tolerance for Earth and Venus, and far 
below for Mars. Manned entry into Jupiter would require a lifting vehicle 
in order to avoid excessive decelerations . 

For vehicles entering from a decaying orbit with aerodynamic lift , 
the maximum heating rate depends strongly on the vehicle weight, shape, 
and dimensions through the parameter w/CpA; maximum heating occurs at a 
velocity of about 0.8 of orbital velocity, and, for any given loading 
W/A, is minimum for entry at C l^ y . This corresponds for common shapes 

to optimum l/d ratios between about 0.5 and 1.0. Because of the 
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coupling between Cp and l/d for any aerodynamic shape , the use of a 
near optimum L/D can reduce the maximum heating rate to no more than 
about one-half that for a nonlifting vehicle. 

The laminar heating rate varies directly as ^W/CpA; hence, by using 
a drag device to increase markedly CpA, such as a drag parachute or flare, 
much larger reductions in heating rate are possible than through the use 
of a trimmed lifting vehicle . 

The total heat absorbed during entry from a decaying orbit increases 
rapidly with lift-drag ratio for vehicles with positive lift. It is a 
minimum for lift-drag ratios near about -0.5, but these negative lifts 
result in excessive decelerations for manned entry into the earth's 
atmosphere; hence the practical optimum for minimizing the total heat 
absorbed in orbital-decay entry of a manned vehicle is near a lift-drag 
ratio of zero. The total laminar heat absorbed, like the laminar heating 
rate, varies directly as ^/W/CpA. 

By inducing entry at a sizable initial angle of descent, the total 
heat absorbed for laminar convection can be reduced substantially. The 
limit of human tolerance to deceleration stress is closely approached for 
nonlifting vehicles entering the earth's atmosphere at an initial descent 
angle of about 3°, under which conditions the total heat absorbed is 0.6 
of that for a decaying orbit having zero initial angle of descent, while 
the decelerations and the maximum heating rate are correspondingly 
increased. However, if a vehicle with small aerodynamic lift (say, 
l/d < 0-7, approximately) enters with a small initial angle, the trajectory 
is a rippling descent which can have a slightly lower maximum heating rate 
as well as smaller total heat absorbed than for gliding entry from a 
decaying orbit. 

The total heat absorbed during the first skip of a lifting vehicle 
entering at a sizable initial angle of descent, is essentially independent 
of both the angle of descent and the velocity of exit from the skip. It 
is a minimum for entry at Ci^y (lift-drag ratios near 0.7). For a given 
w/CpA, this minimum total heat absorbed during the first skip is roughly 
the same as that absorbed during the entry of a nonlifting vehicle enter- 
ing at an initial angle of descent of about 2° . 

In the process of atmosphere braking for stepwise slowing a space 
vehicle from near escape velocity to circular orbital velocity, the 
effects of L/D on peak deceleration, on maximum heating rate, and on 
total heat absorbed are the opposite to the corresponding effects in the 
process of orbital-decay entry. For example, an increase in l/d with a 
given Cp increases the maximum heating rate in atmosphere braking, but 
decreases it in orbital decay. For nonlifting vehicles starting with 
escape velocity and employing atmosphere braking, entry to a planet sur- 
face can be completed on the third pass without the total heat absorbed 
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in any pass exceeding that absorbed for orbital decay, and can be 
completed on the seventh pass without the maximum rate of heating exceed- 
ing that for orbital decay. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. 9> 1958 


16 


NACA TN 4276 


APPENDIX A 

CHECK ON APPROXIMATIONS MADE IN ANALYSIS 


The basic approximation (a) of the analysis, as represented by 
equation (8), can be expressed fairly simply in terms of the transformed 
variable Z and the angle of descent 

[ar/r| _ |u(dy/dt)| = u|sin cp) ^ ^ 

|du/u| |r(du/dt)| ,/pr Z 

Inasmuch as Z/u ~ p , this shows that approximation (a) cannot be valid 
at very high altitudes which are represented by a small neighborhood near 
u = Ui and Z = 0. In figures 28(a) and 28(b), curves of the ratio 
(dr/r)/(du/u) are shown for lifting entry into the earth's atmosphere 
from decaying orbits and for nonlifting entry from deflected orbits with 
various initial angles . It is evident that in the regions near peak 
heating (u = 0.8) and near peak: deceleration (u = 0-3 to 0.5) the basic 
approximation should introduce errors the order of only 1 percent. As a 
vehicle initially enters the atmosphere, however, the decelerations are 
very small and the errors introduced are larger . As a general rule , 
approximation (a) is valid for engineering calculations once the air 
drag has reduced the velocity by about one -half of one percent (see 
appendix B). Approximation (b), that (L/D)|tan cp | « 1 likewise is a 
valid one for heat transfer and deceleration calcinations of vehicles 
with zero or positive lift entering from decaying orbits. As figure 28(c) 
illustrates, approximation (b) may result in substantial errors near 
maximum deceleration for vehicles having negative lift, but still results 
in reasonably small errors near peak heating of such vehicles . 
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APPENDIX B 

MATCHING PRESENT SOLUTION TO KEPLERIAN ELLIPSE 


Let us assume that a retrorocket force , or some other force , has 
deflected an orbiting vehicle into a new Keplerian ellipse which, in the 
absence of drag, would intersect the planet surface at some angle cp G . 

A M zeroth order’ 1 approximation would be to use this angle in the present 
solution as the initial angle cpj_ for the entry. This would be suffi- 
ciently accurate for descent angles greater than a few degrees, but for 
very small angles of descent a more accurate matching of the present 
solution to the Keplerian ellipse may be desirable. 

Since the present solution assumes that 1 dr/r| «| du/u| whereas the 
conservation of angular momentum requires that dr/r = -du/u outside 
the atmosphere, it seems reasonable to select the point of matching where 
the ratio (dr/r)/(du/u) is some value less than unity. Let the descent 
angle at the point of matching be q^, and the velocity be u m . Let us 
confine our attention to a small region near matching, where the density 
is very low, the aerodynamic forces are very small, and the flight path 
is only slightly curved. We represent the Z function in this region 
by the approximation Zj from equation (4o) for constant angle of 
descent, namely, 


Zm = Jpr u m sin q> m ln{xLju. ± ) 

Since u m is only slightly less than ui we approximate In('u m /TTj_) by 
(u m - ui)/ui. Hence from equation (22), it follows that at the matching 
point the ratio r m of terms discarded to terms retained is 

_ dr/r u m sin Cfa „ Uj 

m du/u -/Pr Zm (pr)(uj[ - u m ) 

For Earth pr = 900, so that r m = 1 at u m = 0-999 ^j_, ? m = 0.2 at 
u m = 0.995 ui, and r m = 0.1 at u m = 0-99 ui* Thus, it would be reason- 
able to match the present solution with a Keplerian ellipse at some 
velocity in the range, say, u m = 0-995 Ui to u m = 0.99 The density 

Pm a "k "the matching altitude (from the defining equation (l4) for Z 
and from equation (Bl)) can be determined from 


Pm 

2(m/C D A) 


- Zjn - (u-^ 


- u m) ( - sin cp 0 ) 


(B2) 



^i 

?m 


(-sin cp 0 ) 


Pm - 21 


(B3) 
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or, from p m we can determine the altitude y m and hence the corre- 
sponding value of cpm from the Keplerian ellipse at this altitude • The 
value of cp m so determined would he the value of cpj_ which closely 
matches the present solution for the entry motion. 

An equivalent way of matching would he to select first arbitrarily 
various altitudes y x , y 2 , y 3 , • • • and corresponding densities 
Pl) p 2) p 3) . . . . From the Keplerian ellipse the slightly different 
angles cp 1? cp 2 , cp 3 , . • • could he determined, and hy substitution of 
these into equation (B3) in place of cp 0 , the respective values (m/CpA) 1 , 
(m/Ci>A) 2 , . . . which would bring about proper matching for a given 
value of r m (say, O.l) could be computed. Interpolation would yield 
the matching angle cpm, and the matching altitude y m for any desired 
value of m/CpA. 
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APPENDIX C 

DEVELOPMENT OF SOME APPROXIMATE SOLUTIONS 


The first approximate solution is that for entry of a nonlifting 
vehicle along a spiral path which makes a constant angle cp with respect 
to the local horizontal direction. For this first special case we desig- 
nate the Z function by Zj, and see from the right members of equa- 
tion (19) that 

(sin cp) = 0 = uZj" - s/pr sin cp 
or, after one integration, 

Zj' = -v/pr sin cp In u + constant (Cl) 

The integration constant can be evaluated in terms of the initial veloc- 
ity uq and the angle cp, to yield after one more integration for entry 
from high altitudes (Zi = 0), 


gl 

us/pr sin cp 


u 

= In — =- 


CpApo 


■|3y 


Ui 2mp sin cp 


(C2) 


from which it follows that the dimensionless laminar heating rate 
5 / 2 1/2 

qj = u Zj has a maximum value 


1/4 


(lx ov = 0.247 u ± (pr) n/ sin(-cp) 
■‘•max x 


(C3) 


and the dimensionless total heat absorbed from u = Uj_ to u = 0 is 
evaluated by notin g that the integral of equation (39h) is proportional 
to r(l/2) = Jn/2. 


(Ql) 


Un- 


COS qr 


■Jit/2 


(pr ) 1/4 -s/ sin(-cp) 


(c4) 


This solution for Zj corresponds to setting the left members of 
equation (2l) to zero. In order that the right members of equation (21a) 
also vanish, we see from equation (C2) that this special solution can be 
realized in two ways: (l) by maintaining a true spiral path through pro- 
gramming the lift with velocity in the very special way such that at all 
points 

L (l - u 2 )cos qj 

D u 2 (pr)sin cj> Zn(u/ui) 


(C5) 
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or (2) by entering with a nonlifting vehicle along such a steep path that 
the gravity forces minus the centrifugal forces are negligible compared 
to the vertical component of drag force (this yields essentially a 
straight-line trajectory). Case (l) of spiral trajectories with pro- 
grammed lift, is not easily realized in practice, but case (2) represents 
exactly the physical situation considered by Allen and Eggers (ref. l) 
for their solution to the problem of ballistic entry. Hence, it should 
not be surprising that equation (C2) is identical to their solution. 

This solution for nonlifting vehicles at constant cp does not depart 
significantly from the complete solutions near peak heating (u/uj_ = 0.8) 
except for initial descent angles less than a few degrees, and near maxi- 
mum deceleration (u/ui = 0.4 to 0.6) except for initial descent angles 
less than about 5°* 


As a second special case, we consider smoothly gliding, hypersonic 
flight (u near l) with a large l/d and at sufficiently small descent 
angles that cos cp = 1 and sin cp = cp«L/D. Under these conditions the 
left-hand terms of equation (21a) involving the normal deceleration and 
the vertical component of drag force can be disregarded. The right mem- 
bers yield for the special function Zji representing balance between 
gravity, centrifugal, and lift forces, 


ZH 


1 - u 2 
u N /Pr (L/D) 


(06) 


The flight-path angle is obtained from equation ( 17 ) by differentiating 
equation (c6) 


-cp = 


2 

pru 2 (L/D) 


(C7) 


This particular solution is the same as the solution for gliding flight 
orginally given by Sanger (refs. 2 and 3) for which the aerodynamic heat- 
ing problems have been studied by Eggers, Allen, and Neice (ref. 4). 

This special solution is quite good for L/D ratios greater than about 1 
(for Earth) and hence is adequate for most glide vehicle analysis. It 
cannot be applied, however, to entries with other than zero initial angle, 
inasmuch as extremely small initial angles of descent will result in a 
skipping trajectory for which the vertical acceleration term is not small 
compared to the lift force. For this gliding solution Zjj the maximum 

heating rate proportional to q ma x = (u Z ) max occurs at u = n/2/3 
with 


' 3^(Pr) 1,4 jL7S 


(c8) 


the dimensionless function proportional to the total heat absorbed is 
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Qll = 


(Pr) 


1/4 fz r Ui u 2 du 

>-4> n/i - u 2 



and the range function is 

As II L r Ul u du _ _L_ 1 

r D J o 1 - u 2 2D m l - ui 2 


(CIO) 


as obtained in reference 4. 


As a third special case, we consider entry with lift along a trajec- 
tory wherein the gravity minus centrifugal force is relatively small (see 
eq. (21a)). A skip vehicle, for example, would fall in this category. 

In this case the flight path is determined primarily by a balance between 
the normal acceleration term uZ" , the lift term N /^r(L/D)cos 3 cp, and the 
vertical drag component. The trajectory is, by assumption, influenced 
only secondarily by the gravity minus centrifugal force term 
(l - u 2 )cos 4 cp/(uZ); hence we may render the basic differential equation 
linear by supposing that Z in the denominator of this nonlinear term 
be approximated a priori by some Z function obtained either by neglect- 
ing this nonlinear term or obtained in some other way such as by expand- 
ing Z about uq. By writing cos cp as the "average" value of cos cp 

for the flight path according to the theorem of the mean, we have, after 
one quadrature , 

pU _ 2 

- i = cos 4 cp / 1 ~ ^ du - cos 3 cp -s/pr ^ In = + constant (Cll) 

du u J_ u^Z b ui 

u-i 


at u = u-jj equation (17) shows that (dZ/du) - (z/u) - r sin cpj_ hence 
the descent angle is given by 

JJr (sin cp - sin cpi) = cos 4 cp f — du - cos 3 ^ ^ Jpr In =|~ 

u i 

(C12) 

and the Z function is obtained by solving the first-order differential 
equation (Cll), noting that l/u is an integrating factor. 


Z 

u 


mi A 

— +cos 4 cp 
Ui ' 




(l-u 2 )du , — . , u cos 3 cp 

- — — + «J(3 r sm q>±ln — 


u Z 



u 

Ui 


(C13) 
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By disregarding the gravity minus centrifugal force integral, we obtain 
a special function Zjjj representing balance between normal acceleration, 
vertical drag component, and lift force. 


and 


ZlII 

u 


Zi 

ui 


+ JJr 


u 

sin cp-r Zn — 
1 ui 


cos 3 9 

2 



sin cp = sin cpi - cos 3 cp 



(c lb) 


(C15) 


These last two equations for L/D = 0 reduce to equation (C2). If 
desired, we could substitute Zm (or some other initial estimate of Z) 
into the denominator of the integrand in equation (C13 ) , thus obtaining 
a correction term for the gravity minus centrifugal force term. The suc- 
cess of such a method would depend upon the accuracy and simplicity of 
the initial estimate. 


To illustrate one application of the special solution Zm given 
by equation (Cl4), we consider the first skip only of a lifting vehicle 
entering the atmosphere at a small angle cpi (cos cpi = l) and at orbital, 
velocity (up = l). The first skip is generally the most severe from the 
heat- transfer viewpoint. We have for Z± = 0, 

= n/Pt ^cpiZn u - — Zn 2 u^ (Cl6) 


which can be substituted into the integrand of equation (C13) to yield an 
expression for the gravity minus centrifugal force term. We notice first, 
though, that by definition (Zqn/u) ~ returns to its small initial 
value Zi whenever the vehicle returns to the initial altitude p^. At 
the end of the first skip the velocity is reduced to some value ujjt 
such that 

_ 2c Pi 


in accordance with the results of reference 4. This is the velocity at 
the end of the dip. Since we are considering small angles only, 
-2cp i /(L/D) = 1 - ujjj e , and we may substitute In u = u - 1 in equa- 
tion (Cl6) for the purpose of evaluating the double integral of equa- 
tion (C13) representing the gravity minus centrifugal forces. This 
yields a new Z function 


ziv = (1 - u) 2 + Zm 

u 4^/pr (-cpp ) u 


(C17) 


cpj_2n u - 




= Zi +n/Pt 
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The velocity at the end of the dip is given by 


In u e 


2Ti 

(L/D) + (l/2prcpi) 


(Cl8) 


Since pr = 900 for Earthy the correction term l/(4prcpjj can often be 
disregarded. The path angle is obtained from equation (C15) 


cp = cpi 



(C19) 


so we see that 


cp e 


= cpi - ^ in u e = -cpi 


1 - 


1 + 


1 

2(L/D)prcp i 

1 

2(L/D)prcpi 


(C20) 


If 2(L/D)prcp-[ » 1 these equations reduce to results previously obtained 
by Eggers, Allen, and Neice (ref. 4). In particular, for relatively 
large values of (L/D )cpi , the angle leaving the dip is equal but opposite 
to the initial angle cpi of entry, as noted in reference 4. After a 
skip, a period of weightlessness follows at an essentially constant 
velocity u e under conditions where the vertical acceleration 
(dv/dt) = g - (u e 2 /r) is constant; hence the duration of weightlessness 
At = 2v e /(dv/dt) is 


At 


2ug (cpg ) 

g(l - u e 2 ) 


(C21) 


After this period, a second entry occurs at nearly the same angle as the 
first entry, only at the reduced velocity u e . The maximum laminar- 
heating rate occurs near, but not at the bottom of the dip (cp = 0) at 
which point the velocity u m is given by In u m = cpi/(L/D). By substi- 
tuting this into equation (Cl6) the approximate maximum laminar heating 
rate is then represented by 

acpi 


5/2 1/2 1/4 l/d ( _< Pi) 

= (u m Zfa ) = (pr) e 


(C22) 


, N l/4 / \ 

(pr) (-9i) 


•J 2L/D 


for -cpi « p 


An interesting result concerns the total heat absorbed in a single 
skip starting from satellite velocity (up = l)» The total heat absorbed 
is obtained from equation (39) together with 
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Q = 


, T 3/2 

1 u du 


42 


u du 


u e 


4z 


(er)1, ‘^' ue 


By employing the same approximation Jn u = -1 + u the integral can be 
evaluated . 


Q = 


42 

(Pr) 1/4 JI^ 


it 1 + 


l/d 


(C23) 


This is essentially independent of cpi, since cpi/(L/D) for many skip 
vehicles would be small. 


Q = 


42 * 

(pr) 1 Vl/d 


for -cp^ 



Although the maximum heating rate in a skip is proportional to the initial 
angle of descent cp^ the total heat absorbed is essentially independent 
of both the_initial angle cp-p and the exit velocity u e of the skip. 
Since Q ~ Q,/ S /Cp ; we see that Q ~ 1 /JCl, which means that the least heat 
is absorbed by skipping at CL max - For inclined flat surfaces in hyper- 
sonic flow, simple Newtonian theory yields CL max = 0-77 at an angle of 
attack of 55° > at which angle L/D = 0.71* Hence 

(Q) mi n Q = 42 ,/( P r) 1/4 ^7l = 0.96 


This value is compared elsewhere in the DISCUSSION with corresponding 
values of Q for other types of entry. 


If a skip vehicle does not enter initially at orbital velocity, but 
at some different value uq, then the corresponding equations with gravity 
and centrifugal forces neglected indicate the bottom of the dip to be at 
a velocity u m given by Zn(%]/ui) = cp-j_/ (L/D ) . At this point the heating 
rate is represented by 


dm = 


acpi 

x 1 / 4 - 3 L / D ( -c Pi ) 
(|3r) u ± e 

^2L/D 


(C24) 


The exit of the skip occurs at a velocity u e given by 

Zn(u e /ui) = 2cpi/ (L/D) 


The dimensionless total heat absorbed is approximately 


42 ira± 2 

(Pr) 1/4 «yL/D 


(C25) 
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APPENDIX D 

INTEGRATION OF BASIC NONLINEAR EQUATION 


Many numerical methods could be used to compute stepwise a Z func- 
tion from a nonlinear equation such as 

*7 ~\ T 

uZ ,! - Z T + — = — ' cos 4 cp - — cos 3 cp (Dla) 

where Z r - ~ = V^r* sin cp. A study has not "been made of the best way to 
u 

integrate such an equation, or of whether or not an alternate form of this 
equation, such as 



f (=2 - 1 ) cos 4 cp - ^ cos 3 cp 


(Dlb) 


where F = Z/u and uF T = 0r sin cp, may he preferable for purposes of 
integration. The particular method employed, while probably not extremely 
accurate, is simple in the sense that it involves merely the repetition 
of a large number of identical operations. Suppose we know at some ini- 
tial point u n the values of and Z n T . Then from the differential 
equation we have for the second derivative (with cos cp set equal to 
unity for purposes of simplicity in illustrating the method). 


V =n*(V 


and for the third derivative 


— + ~ u n 2 

u n ~n^n 



(D2) 


Zn\ _ 1 - Un 2 Z n ' 1 + un 2 
Un J u n Zn 2 u n 2 Z n . 

Hence a Taylor expansion for Zn +1 and Zn +1 ' at the next point 
yields 


Zn 


Tit _ 


u n 


1 ( Z n ' 
_u n 


(D3) 

u n +i 


Zn+i 

Zn+i' 


Zn + Z n ' Ai + Zn" 


(Au) S 

2 


+ Zn"' 


(An) 3 

6 




... (ast 


+ Zn Au + Zn 


2 
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while the above equations (D2) and (D3) yield Z n +i" and Zn+i m when n 
is replaced by n+i in the formulas. Thus the process can be continued. 
For most cases the Z functions are fairly smooth, and the inclusion ^ 
of Z"' is unnecessary in the above procedure if sufficiently small Au 
are used. For the present calculations Z m was omitted; Ad = 0.001 was 
employed for u^ = 1, and Au = 0.002 for ui = 1.4 and uj[ = 1.2. For skip 
vehicles, the Z function can vary quite rapidly and^the inclusion 
of Z'" presumably would enable larger increments Ad to be used. 

This particular procedure requires a knowledge of nonzero values 
Z o and Zq' at some initial point Uq. Hence the first step is taken 
analytically. For decaying orbits an analytical representation in the 
vicinity of u = 1, where (l - u 2 )/u = 2(l - u), is 


- 2 J (1 - 5o f /2 

(D4) 

= -46 (1 - Uo) 1/2 

(D5) 


since these yield Z 0 " = 2(l - u 0 )/Z 0 and correspond to values of both 
Z 0 ' and Z 0 /uo small compared to Z 0 " (see eq. (Dl)). Equations (D4) 
and (D5) would apply to a lifting vehicle provided (l - u G ) is selected 

T , 1 y p 

small enough so that n/ (3r — « (l - u G ) . It the L/D ratio is 

D 

large, we can use the Zjj function to obtain 


Zo 


z 


T 

o 


1 - Up 2 
\t~$r (L/D)u 0 

1 — 2 

1 + Up 

■v/F (L/D)u 0 2 


(D6) 

(D7) 


For re-entry with an initial angle cp^ at initial velocity u^ we can 
use the Zjjj function for the first step. 


Z 0 = n/F Uo 


U 


sin cp 4 In — - 
± u*? 


cos^cpi ( L\ ln 2 _u 


p. 


u i. 


V 


= F 


sin au + — 
1 u~ 


(D8) 

(D9) 
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TABLE I.- VALUES OF Z FUNCTIONS AND RELATED QUANTITIES FOR 


L/D = 0 AND vl ± = 1 



Z 

_Cp Earth 

deg 

Q 

(-) 

\ r 'Earth 

t Earth 

sec 

Z 

“ cp Earth 

deg 

$ 

/Earth 

^ Earth 
sec 

z 

_C ^Earth 

deg 

Q 

(*) 

\ r /Earth 

t Earth 

sec 

u 



= 0° 

irth 


Earth = 

0.5° 


(-9. 

^ Earth 

: 1° 


0.995 

0.00058 

0.33 

0 

0 

0 

0.00131 

0.50 

0 

0 

0 

0.0026 

1.00 

0 

0 

0 

.99 

.00165 

.48 

.156 

.169 

140 

.00270 

.57 

.113 

.088 

72 

.0053 

1.04 

.080 

.044 

37 

• 98 

.00467 

.68 

• 338 

.288 

238 

.00603 

.71 

.265 

.169 

140 

.0108 

i.n 

.191 

.088 

73 

.96 

.01315 

.96 

.550 

.372 

309 

.01457 

.95 

.460 

.240 

199 

.0229 

1.26 

.342 

.129 

108 

.94 

.0241 

1.19 

.688 

.409 

341 

.0253 

1.16 

.593 

.274 

229 

.0361 

1.40 

.451 

.152 

128 

.92 

.0369 

1.39 

•792 

.431 

360 

.0378 

1.36 

.695 

.296 

248 

.0505 

1.53 

• 537 

.168 

141 

.90 

.0515 

1-57 

.875 

.446 

374 

.0519 

1.53 

.777 

.311 

262 

.0660 

1.67 

.609 

.180 

152 

.85 

.0939 

1.98 

1.030 

.470 

396 

.0934 

1.94 

.932 

.335 

284 

.1092 

2.01 

.750 

.199 

170 

.8 

.1435 

2.36 

1.140 

.485 

411 

.1421 

2.33 

1.042 

.349 

298 

.1580 

2.34 

.853 

.212 

183 

• 75 

.1991 

2.73 

1.223 

.494 

421 

.1970 

2.70 

1.126 

.359 

309 

.212 

2.69 

.933 

.221 

192 

.7 

.260 

3.H 

1.288 

.502 

429 

.257 

3.08 

1.191 

.366 

317 

.270 

3.05 

.996 

.228 

200 

.65 

.324 

3.50 

1.340 

.507 

436 

.321 

3.48 

1.243 

.372 

324 

.333 

3.42 

1.047 

.233 

207 

.6 

.392 

3.91 

1.381 

• 512 

442 

.389 

3.89 

1.285 

.377 

330 

.398 

3.83 

1.088 

.238 

213 

• 55 

.463 

4.36 

1.415 

.516 

448 

.460 

4.35 

1.318 

.381 

336 

.466 

4.28 

1.122 

.242 

218 

.5 

.536 

4.86 

1.442 

.519 

453 

.533 

4.85 

1.346 

.384 

341 

.537 

4.77 

1.149 

.245 

223 


.610 

5.43 

1.464 

.522 

458 

.607 

5.42 

1.367 

.387 

346 

.608 

5.34 

1.170 

.248 

228 

.4 

.684 

6.09 

1.481 

.525 

463 

.681 

6.08 

1.385 

.390 

351 

.680 

6.00 

1.I88 

.251 

233 

• 35 

• 757 

6.89 

1.495 

.527 

468 

.754 

6.89 

1.399 

.392 

356 

• 752 

6.81 

1.202 

.253 

238 

.3 

.827 

7-89 

1.505 

.529 

473 

.824 

7-89 

1.409 

.394 

361 

.821 

7-82 

1.212 

.255 

244 

.25 

.892 

9.21 

1.513 

.531 

479 

.890 

9.22 

1.417 

.396 

367 

.885 

9.15 

1.220 

.257 

249 

.2 

.949 

11.09 

1.519 

.533 

486 

.947 

11.10 

1.423 

.398 

374 

.942 

n.o4 

1.226 

.259 

256 

.15 

.992 

14.04 

1.523 

.535 

493 

.991 

14.05 

1.427 

.399 

382 

.986 

14.01 

1.230 

.260 

264 

.1 

1.009 

19.52 

1.525 

.536 

504 

1.009 

19.53 

1.429 

.401 

392 

1.005 

19.52 

1.232 

.262 

275 

.05 

• 958 

33.3 

1.527 

.538 

521 

.958 

33.3 

1.431 

.402 

409 

.956 

33.3 

1.234 

.263 

292 

.025 

.825 

50.7 

1.527 

• 538 

537 

.825 

50.7 

1.431 

.403 

425 

.824 

50.8 

1.234 

.264 

307 



(-q>i) E 

arth = 2< 

3 


Earth 

= 3° 


Earth 1 

= 4° 


0.995 

.0052 

2.00 

0 

0 

0 

.0078 

3.00 

0 

0 

0 

.0105 

4.00 

0 

0 

0 

• 99 

.0105 

2.02 

.057 

.0222 

IB 

.0157 

3.02 

.047 

.0148 

12 

.0209 

4.01 

.040 

,om 

9 

.98 

.0210 

2.06 

.136 

.0443 

37 

.0313 

3.04 

.111 

.0295 

24 

.0416 

4.03 

.097 

.0222 

18 

.96 

.0422 

2.13 

.246 

.0663 

55 

.0624 

3.09 

.201 

.0443 

37 

.0827 

4.07 

.175 

.0333 

28 

.94 

.0638 

2.21 

.327 

.0790 

66 

.0934 

3.15 

.268 

.0529 

44 

.1233 

4.11 

.233 

.0398 

33 

• 92 

.0858 

2.29 

.393 

.0880 

74 

.1243 

3.20 

.323 

.0591 

50 

.1633 

4.15 

.281 

.0445 

38 

.90 

.1080 

2.37 

.449 

.0949 

80 

.1550 

3.26 

.370 

.0639 

54 

.203 

4.19 

.322 

.0481 

41 

.85 

.1651 

2.59 

.560 

.1073 

92 

.231 

3.41 

.463 

.0726 

62 

.300 

4.31 

.404 

.0548 

47 

.8 

.224 

2.82 

.646 

.1159 

100 

.306 

3.57 

.536 

.0788 

68 

• 393 

4.44 

.469 

.0596 

52 

• 75 

.284 

3.07 

.714 

.1225 

107 

.380 

3.76 

.595 

.0837 

74 

.482 

4.58 

.521 

.0634 

56 

.7 

.346 

3.34 

.769 

.1278 

113 

.452 

3.96 

.643 

.0877 

78 

.568 

4.74 

.564 

.0666 

59 

.65 

.409 

3.64 

.815 

.1322 

119 

.523 

4.20 

.683 

.0911 

82 

.649 

4.93 

.599 

.0693 

63 

.6 

.474 

3-97 

.852 

.1360 

124 

.591 

4.46 

.717 

.0941 

86 

.726 

5.14 

.630 

.0718 

66 

• 55 

• 538 

4.35 

.883 

.1393 

128 

.657 

4.77 

.744 

.0968 

90 

• 797 

5.39 

.655 

.0739 

69 

.5 

.603 

4.79 

.908 

.1422 

133 

.721 

5.13 

.767 

.0992 

94 

.863 

5.69 

.676 

.0759 

72 

.45 

.668 

5.30 

.929 

.1448 

137 

.781 

5.56 

.787 

.1014 

98 

.922 

6.06 

.693 

.0778 

75 

.4 

.732 

5.91 

.946 

.1472 

142 

.838 

6.09 

.802 

.1035 

102 

.974 

6.51 

.708 

.0795 

79 

• 35 

.794 

6.67 

.959 

.1493 

147 

.890 

6.76 

.815 

.1054 

106 

1.018 

7.08 

.719 

.0812 

82 

.3 

.853 

7.63 

.969 

.1513 

152 

.936 

7.63 

.824 

.1072 

no 

I.052 

7-85 

.729 

.0828 

86 

.25 

.907 

8.93 

• 977 

.1532 

157 

.975 

8.82 

.832 

.1089 

115 

1.076 

8.92 

.736 

.0843 

91 

.2 

.954 

10.79 

.983 

.1550 

164 

1.005 

10.57 

.838 

.1106 

121 

1.086 

10.52 

.741 

.0859 

97 

.15 

.989 

13.75 

.987 

.1567 

172 

1.021 

13.42 

.841 

.1122 

129 

1.079 

13.20 

.745 

.0874 

104 

.1 

1.002 

19.28 

.989 

.1583 

182 

1.016 

18.87 

.844 

.1137 

139 

1.048 

18.45 

.747 

.0889 

n4 

.05 

.952 

33.2 

.991 

.1598 

200 

.952 

32.8 

.845 

.1153 

157 

.960 

32.3 

.748 

.0904 

131 

.025 

.822 

50.8 

• 991 

.1605 

215 

.820 

50.6 

.845 

.1160 

172 

.821 

50.3 

.74? 

. 09 n 

146 


NACA TN 4276 






6o 


NACA TN ^276 





_ 





Altitude f t x 10' 



o 


Figure 1.- Comparison of exponential approximation with ARDC model of Earth atmosphere ( 1956 ). 
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Dimensionless parameter 
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Figure 2.- Dimensionless parameter 

atmosphere. 


for ARDC model of Earth 
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Trajectory 

Portion of trajectory over 
which analysis is applicable 







Entry from deflected 
orbit with initial angle (p f 
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Figure 3 .- Sketches of typical entry trajectories and. portions to which 

present analysis applies . 
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(a) Nonlifting vehicles. 

Figure 4.- Values of Z functions for entry from decaying orbits into 

planetary atmosphere. 
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(b) Vehicles with L/D > 0. 


Figure 4.- Concluded. 





Ch 



Figure 5.- Comparison of Z functions for nonlifting entry with analysis of Allen and Eggers; 

u ± =0.9 (23,400 fps for earth). 
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Figure 


30uZ 



(a) Nonlifting vehicles. 


6. - Values of Z functions for entry from orbital velocity at 
Initial angles of descent; = 1. 


68 


NACA TN 14-276 



f’igure 6.- Continued. 
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Dimensionless velocity, u 

(c) Vehicles with L/D = 0.5. 


Figure 6.- Continued. 
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Dimensionless velocity, D 


(d) Vehicles with L/D =0.7 
Figure 6.- Continued. 
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Dimensionless velocity, u 

(e) Vehicles with L/D = 1. 


Figure 6.- Concluded. 
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vehicles; Uj_ = 1.4. 


Figure 7. - Values of 


Angle of descent, -(p,deg 
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0 .2 .4 .6 .8 1.0 

Dimensionless velocity, u=jj~ 

(a) Angle of descent and altitude -velocity trajectory. 

Figure 8.- Comparison of present approximate analysis with more exact 
machine calculations for ARDC model atmosphere. 
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(t>) Deceleration and distance traveled. 
Figure 8.- Concluded. 
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rresem analysis, using vpr 

Present analysis, using semi-local value 

of «/JTr for ARDC atmosphere (fig. 2 ) 



• Numerical integration for ARDC atmosphere 

( Rubesin-Goodwin) . 

1 1 1 L_ l 1 1 1 
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Figure 9.- Example application of present analysis to a model atmosphere, 
and comparison with more exact machine calculations; nonlifting entry 
from decaying orbit. 


Horizontal deceleration, gs 
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Figure 10 


Effect of lift -drag ratio on deceleration for entry into 
Earth atmosphere from decaying orbits. 


Maximum deceleration in earth g's 
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Lift-drag ratio, 


Figure 11.- Effect of lift-drag ratio on maximum deceleration for entry 
into various planetary atmospheres from decaying orbits. 


Maximum resultant deceleration, earth g's 
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Figure 12.- Comparison of decelerations and duration for entry into 
various planetary atmospheres from decaying orbits. 
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Figure 13. - Reynolds number at peak heating for entry from decaying orbits 

into Earth atmosphere. 
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Figure 14.- Effect of lift-drag ratio on laminar heating rate for entry 

from decaying orbits. 


Maximum heating rate 
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Lift-drag ratio (b) 

Earth 


Figure 15.- Effect of lift-drag ratio on maximum laminar heating rate at 
stagnation point for entry from decaying orbits. 


for earth, °R ft 



Figure 16.- Maximum radiation -equilibrium, temperature at laminar stagnation point for entry from 

decaying orbits into Earth atmosphere. 
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Figure 17.- Maximum surface temperature for entry into various planets 
from decaying orbits (u-^ =1, cp-p = 0). 
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Figure 18.- Effect of lift -drag ratio on total heat absorbed during entry 
from decaying orbits into Earth atmosphere. 
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Lift-drag ratio, ^ 

Figure 19.- Total heat absorbed during entry into various planets from 

decaying orbits (u^ = 1, cp^ = 0). 


Maximum deceleration, g's 
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Figure 20.- Effect of initial angle of entry on deceleration, laminar 
aerodynamic heating, and range of nonlifting vehicles. 


Dimensionless total heat absorbed, 
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(b) Maximum laminar heating rate and total heat absorbed. 


Figure 20.- Continued 


Distance, statute miles 
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Initial angle, - tp j, deg 


(c) Horizontal distance and horizontal miss distance. 
Figure 20.- Concluded. 
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Figure 21.- Reynolds number at peak heating for nonlifting entry from 
deflected orbits into Earth atmosphere. 


Dimensionless moximum heating rate, 
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(a) Maximum laminar heating rate. 

Figure 22.- Effect of initial angle of entry on laminar aerodynamic 
heating and range of lifting vehicles. 
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Figure 22. - Continued 




Initial angle, -((p|) Earth ,deg 

(c) Entry range from u = 0.99 to u = 0. 


Figure 22. - Concluded 


Dimensionless total heat absorbed during first skip, 
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Figure 23.- Dimensionless maximum heating rate and total heat absorbed 
for a composite entry with L/D = 0 down to velocity uv,, and then 
L/D =0.7 thereafter. 
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(a) Maximum horizontal deceleration. 

Figure 2k. - Atmosphere braking for single pass starting with Uj_ = 1.4. 
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(t>) Maximum heating rate 
Figure 24. - Continued. 


Dimensionless heat absorbed in single pass, 
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(c) Heat absorbed in single pass. 
Figure 24.- Concluded. 
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(a) Maximum horizontal deceleration and maximum heating rate. 

Figure 25.- Atmosphere braking for single pass starting with uj; = 1.2. 
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Dimensionless velocity at exit of braking, u 

(b) Heat absorbed in single pass. 
Figure 25. - Concluded. 
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Figure 26.- Maximum temperature during single atmosphere braking pass 

starting with Uj_ = 1.4. 
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Figure 27. - Heat absorbed during single atmosphere braking pass starting 

with Uj_ = 1.4. 
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(a) Approximation (a) for various L/D. 
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Figure 28.- Check on approximations made in analysis. 
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